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Abstract 
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. Context. Molecular counterparts to atomic jets have recently been detected within 1000 AU of young stars at early evolutionary 

' stages. Reproducing these counterparts is an important new challenge for proposed ejection models. 

' Aims. We explore whether molecules may survive in the magneto-hydrodynamic (MHD) disk wind solution currently invoked to 

reproduce the kinematics and tentative rotation signatures of atomic jets in T Tauri stars. 

Methods. The coupled ionization, chemical, and thermal evolution along dusty flow streamlines is computed for the prescribed MHD 
disk wind solution, using a method developed for magnetized shocks in the interstellar medium. Irradiation by (wind-attenuated) 
coronal X-rays and far-ultraviolet photons from accretion hot spots is included, with an approximate self-shielding of H2 and CO. 
Disk accretion rates of 5x 10"^, 10"* and lO^^Moyr"' are considered, representative of low-mass young protostars (so-called 'Class 0'), 
. evolved protostars ('Class F) and very active T Tauri stars ('Class IF) respectively. 

Results. The disk wind has an "onion-like" thermo-chemical structure, with streamlines launched from larger radii having lower 
temperature and ionisation, and higher H2 abundance. The coupling between charged and neutral fluids is sufficient to eject molecules 
from the disk out to at least 9 AU. The launch radius beyond which most H2 survives moves outward with evolutionary stage, from 
^ 0.2 AU (sublimation radius) in the Class disk wind, to 1 AU in the Class I, and > 1 AU in the Class II. In this molecular 
^ ■ wind region, CO survives in the Class but is significantly photodissociated in the Class I/II. Balance between ambipolar heating and 

molecular cooling establishes a moderate asymptotic temperature ^ 700 - 3000 K, with cooler jets at earlier protostellar stages. As a 
result, endothermic formation of H2O is efficient, with abundances up to ^ lO""*, while CH^ and SH^ can reach > lO"*" in the hotter 
and more ionised Class I/II winds. 

^ ■ Conclusions. A centrifugal MHD disk wind launched from beyond 0.2 - 1 AU can produce molecular jets/winds up to speeds 

QQ ' 100 km s"' in young low-mass stars ranging from Class to active Class II. The models predicts a high ratio of H2 to CO and an 

I increase of molecular launch radius, temperature, and flow width as the source evolves, in promising agreement with current observed 

, trends. Calculations of synthetic maps and line profiles in H2, CO and H2O will allow detailed tests of the model against observations. 

, Key words, astrochemistry - stars: formation - stars: mass loss - ISM: jets and outflows - ISM: molecules - Accretion, accretion disks 
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1 . Introduction all three processes probably coexisting to some degree in young 

stars (see Ferreira et al. 2006 and refs. therein). 

Powerful supersonic bipolar jets are ubiquitous in young accret- Qf the wide variety of models pubUshed, stationary solutions 

ing stars, suggesting that jets may play a key role in extract- for "X-winds" and extended MHD disk-winds are of particular 

ing excess angular momentum from the accretion flow. The high interest as they provide synthetic predictions in atomic lines that 

^ . ejection efficiency, and the universal jet collimation properties reproduce several aspects of jet colUmation and acceleration in 

independendy of circumstellar envelope density, both strongly microjets from T Tauri stars (Shang et al. 1998; Shang et al. 

favor a magneto-hydrodynamical (MHD) ejection and collima- 2002; Cabrit et al. 1999; Garcia et al. 2001b; Pesenti et al. 2004; 

tion process (see e.g. Konigl & Pudritz 2000; Shu et al. 2000; pyo et al. 2003, 2006; Cabrit 2007b). Extended MHD disk winds 

Cabrit 2007a). However, the main jet origin is still the subject have received further attention since the tentative detections of 

of intense debate. Proposed options include a self-collimated, jgt rotation signatures, that suggest an outer launch radius for the 

steady centrifugally-driven disk wind — either from the inner atomic jet' in T Tauri stars ^0.2-3 AU, well beyond the disk 
disk edge near corotation ("X-wind") or from a more extended 

range of radii — a pressure-driven stellar wind, and unsteady ' We will use the term "atomic jet" rather than "optical jet" to refer 

ejections by reconnexions in the sheared stellar magnetosphere, to the (mosfly neutral) fast jet component emitting in atomic/ionic lines, 

as infrared detections e.g. in [Feii] are increasingly frequent (see Bally 

et al. 2007, for a review). Atomic jets differ by their high speed and 

Send offprint requests to: panoglou@phl . uni -koeln . de colUmation from the "neutral atomic outflows" detected in H i in a few 
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* Present address: I. Physikalisches Institut, Universitat zu Koln, young protostars, whose bulk mass is moving at < 20 km s (Lizano 
Zulpicher StraBe 77, 50937 Koln, Germany et al. 1988; Rodriguez et al. 1990). 
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inner edge (Bacciotti et al. 2002; Anderson et al. 2003; Pesenti 
et al. 2004; Coffey et al. 2004, 2007). Evidence for depletion of 
refractory species in atomic jets also suggests ejection from rela- 
tively large disk radii, outside the dust sublimation radius (Nisini 
et al. 2005; Podio et al. 2006, 2009; Agra-Amboage et al. 201 1). 
If they remove most of the angular momentum from the accre- 
tion flow, such extended MHD disk winds would have a strong 
impact on the conditions for planet formation and migration in 
disks (Combet & Ferreira 2008). Therefore it is important to find 
further observational tests of their presence and radial extent. 

An important new test for proposed jet engine(s) is whether 
they can explain the molecular jet/wind counterparts recently de- 
tected within 1000 AU of young stars. Molecular jets were first 
discovered emanating from the so-called 'Class 0' protostars at 
the earliest stage of star formation, in CO, SiO, H2, and H2O 
masers (Bachiller & Cernicharo 1990; Guilloteau et al. 1992; 
Gueth & Guilloteau 1999; McCaughrean et al. 1994; Claussen 
et al. 1998). These narrow molecular jets are clearly distinct 
from the "standard" swept-up molecular outflows in that they 
are much faster, with deprojected speeds of 60 - 150 km s"' 
as opposed to < 20 km s"' in the swept-up gas, much more 
highly collimated on-axis, and with specific chemical enhance- 
ments (Bachiller & Cernicharo 1990; Gueth & Guilloteau 1999; 
Santiago-Garci'a et al. 2009; Tafafla et al. 2010). Their colU- 
mation scale, ejection/accretion ratio, and variation of physical 
parameters with velocity are strikingly similar to those of atomic 
T Tauri jets, pointing to a similar origin (Cabrit et al. 2007; Lee 
et al. 2007; Nisini et al. 2007). Counterparts in H2, and more 
rarely in CO, were then identified towards more evolved infrared 
protostars of Class I (Cernicharo & Reipurth 1996; Davis et al. 
2001b, 2002; Davis et al. 2006; Bally et al. 2007) and optically 
visible T Tauri stars of Class II (Takami et al. 2004; Takami et al. 
2007; Herczeg et al. 2006; Pety et al. 2006; Beck et al. 2008). 
They consist of high-speed H2 in the jet beam itself, and/or of 
slow wind/cavities at ^ 30 km s"' encompassing the fast ionic 
jet, on scales of a few 100 AU. The molecular emission tends to 
peak at lower velocities than atomic lines (Davis et al. 2003). 

Three processes are generally invoked to produce molecular 
jets/winds in young stars: (i) entrainment of ambient molecular 
material, (ii) gas-phase reformation of H2 in a dust-free atomic 
wind (e.g. stellar or magnetospheric); (iii) ejection/reformation 
of molecules in a dusty wind from the disk surface. 

The entrainment hypothesis is intuitively appealing, but ap- 
pears quite challenging. Turbulent mixing of ambient gas along 
the walls of a supersonic jet beam predicts too low columns 
of H2 (Taylor & Raga 1995). Shock entrainment in narrow 
bowshock wings along the jet is potentially more efficient, but 
would require that ambient molecular gas can quickly mix- 
in and refill the cocoon behind preceding bowshocks (Davis 
et al. 2001a; Raga & Cabrit 1993). Swept-up molecular cav- 
ities driven into ambient gas by an atomic wide-angle wind 
have also been proposed in Class II sources where molecular 
counterparts are broader than the atomic jet (see e.g. Takami 
et al. 2007). However, the nature of such a wide-angle com- 
ponent remains unclear; the fast ^ 100 km s ' and powerful 
wide-angle wind predicted by the "X-wind" model, propagating 
into a flattened envelope, would sweep up a cavity of width to 
length > 1 at 200 AU from the source in less than 1000 years 
(Lee et al. 2001; Shang et al. 2006). This is difficult to rec- 
oncile with the moderate width/length < 1 of molecular "cav- 
ities" in Class II sources of age ^ 10^ yrs, unless the wide- 
angle wind is much slower/weaker (e.g. a disk wind launched 
from large radii). In any case, an ambient swept-up cavity would 
not easily explain the narrow, high- velocity molecular jets seen 



in Class O/I sources. The presence of blueshifted H2 or CO up 
to 50 - 100 km s'^ within 10 AU of Class II / Class I sources 
(Herczeg et al. 2006 and Herczeg, private communication) also 
seems difficult to explain with entrainment. 

Alternatively, it has been argued that molecular counter- 
parts may trace the primary ejected material itself (Cernicharo 
& Reipurth 1996; Takami et al. 2004; Pety et al. 2006; Hirano 
et al. 2006). Studying the formation/survival of H2 in protostel- 
lar winds is a difficult problem as it requires the combination of 
state of the art MHD wind models with chemical networks and 
realistic radiation fields and heating/cooling terms, few authors 
having addressed it self-consistently in the past. 

Reformation of H2 in dust-free stellar winds by ion reactions 
has been the most extensively studied so far. Spherical models 
with prescribed temperature laws showed that H2 forms and sur- 
vives only at high mass-flux rate and if temperature drops rapidly 
below 2000 K (Rawlings et al. 1988; Glassgold et al. 1989; 
Rawlings et al. 1993). Ruden et al. (1990) solved in parallel for 
the ionisation, H2 abundance, and thermal evolution, including 
ion-neutral "drag" heating. They found that less dense winds are 
hotter and more collisionally dissociated, so that the final H2 
abundance exceeds a few percents only at very high wind mass- 
fluxes > 10"''Mo yr"'. Such a model does not seem able to ex- 
plain the presence of detectable H2 in low-luminosity Class O/I/II 
jet sources, where total mass-flux rates are 1-3 orders of mag- 
nitude smaller (Lee et al. 2007; Antoniucci et al. 2008; Hartigan 
et al. 1995). Far-ultraviolet excess from accretion further reduces 
the predicted H2 (and SiO) abundance at low wind mass-flux, 
even though CO is efficiently formed (Glassgold et al. 1991). 
Alternatively, Raga et al. (2005) show that Ht abundances of up 
to 10% could be reached behind dense internal shocks in a dust- 
free jet. However, the jet velocity variability amplitude should be 
less than 15 km s"', which is quite restrictive. The effect of mag- 
netic cushioning and of far-ultraviolet (FUV) photons remains to 
be investigated. 

Here we investigate the third possibility, namely whether 
molecular jet counterparts could be tracing dusty disk winds. 
We focus on self-similar, centrifugally-driven MHD disk winds. 
Pioneering calculations by Safier (1993) indicated that H2 could 
survive collisional dissociation in the outer regions of such disk 
winds, beyond a minimum launch radius ~ 0.5 - 3 AU depend- 
ing on the adopted MHD solution. However, the rate of colli- 
sional dissociation may have been underestimated, as the ion- 
isation fraction was obtained assuming atomic gas. The faster 
dissociative recombinations occuring in molecular gas would en- 
hance ion-neutral drag, increasing the wind temperature and the 
dissociation rate. The calculations also did not consider H2 de- 
struction by stellar FUV photons, coronal X-rays, or endother- 
mic reactions with O and OH. 

In the present paper, we explore this problem one step further 
by integrating in parallel for the detailed chemistry, ionization, 
and thermal state along MHD disk wind streamlines. We use a 
model developed for magnetized molecular interstellar shocks 
with an extended chemical network of 134 species (Flower & 
Pineau des Forets 2003), adding irradiation by (attenuated) coro- 
nal X-ray photons and FUV photons from the accretion shock. 
Another difference with respect to Safier (1993) is the use of 
a slower and denser MHD accretion-ejection solution that better 
matches the rotation signatures, poloidal velocities, and ejection- 
accretion ratio in atomic T Tauri jets (Pesenti et al. 2004; Ferreira 
et al. 2006). The flow dynamics, thermo-chemistry, and treat- 
ment of FUV and X-ray radiation, are described in detail in 
Section 2, which may be skipped in a first reading. Our results 
are presented in Section 3 for a range of accretion rates and stel- 
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lar masses representative of the Class 0, Class I, and most ac- 
tive Class II phases of solar-type stars. In Section 4, the results 
are compared to static irradiated disk models and to observa- 
tional trends, and model limitations are discussed. Conclusions 
are summarized in Section 5. 



2. Model description 

2.1. Dynamical MHD solution 

The adopted dynamical MHD model sets the total density, veloc- 
ity, magnetic field and current in space. A large scale magnetic 
field of bipolar topology is assumed to thread the disk on a large 
radial extension. The wind model satisfies the full single-fluid 
MHD equations and has been obtained by assuming steady-state, 
axisymmetry, and a self-similar functional form for all physical 
variables. These models are an extension of the Blandford & 
Payne (1982) ideal MHD jet models in that they also solve for 
the full 3D dynamics of the underlying resistive keplerian accre- 
tion disk (where accreting matter needs to cross field lines), in- 
cluding all dynamical terms in the vertical, azimuthal, and radial 
directions. In particular, the disk compression and (slight) sub- 
keplerian rotation induced by the magnetic field, addressed in 
Shu et al. (2008), are accurately treated. The accretion-ejection 
transition, wind mass-loading, and slow point crossing, are thus 
all determined self-consistently, given a vertical profile of tem- 
perature and resistivity in the disk (see Ferreira 1997; Casse & 
Ferreira 2000 for details). 

A key free parameter of the models that is most amenable to 
observational constraints is ^, the disk ejection efficiency, fixing 
the accretion rate as a function of the radius as Mace ■ The to- 
tal mass outflow to disk accretion rate over a radial extension ri„ 
to Tout is then IMj/M^cc = ^ln(rout/'"in)- The parameter ^ is also 
linked to the jet physics. As long as wind magnetic torques dom- 
inate, it is inversely related to the Blandford & Payne magnetic 
lever arm parameter A through; 



A = 



1 + 



1 

2^' 



(1) 



where ro is the midplane radius of the magnetic surface along 
which the matter flows and is its cylindrical radius at the 
poloidal Alfven point (Ferreira 1997). The maximum asymp- 
totic jet velocity and specific angular momentum along each 
streamline are given respectively by ^ yfGMZJro V2/1 - 3 
and n'^ ^ AyjCM^^rQ (Blandford & Payne 1982), and are thus 
smaller for decreasing A and increasing ejection efficiency ^. 

In Garcia el al. (2001a), we calculated the thermal-ionisation 
structure for atomic MHD disk wind solutions with high lever 
arms (42 < A < 84), obtained for a vertically isothermal disk 
profile (Ferreira 1997). Garcia et al. (2001b) presented synthetic 
forbidden line predictions and showed that the predicted line 
centroid velocities were on average too high when compared to 
observations, and that models with lower lever arms should be 
favored. This conclusion was later confirmed by the small rota- 
tion signatures observed in T Tauri atomic jets, which rule out 
values of A > 20 (Pesenti et al. 2004; Ferreira et al. 2006). 

Therefore, in the present study we use a "slow" disk wind 
solution with a smaller lever arm parameter A - 13.8, developed 
by Casse & Ferreira (2000). This particular value was chosen as 
best reproducing the possible rotation signatures reported across 
the outer layers of the DO Tau atomic jet, with an inferred outer 
launch radius rout - 3 AU (see Pesenti et al. 2004). The same 
model can also account for the highest velocity ^ 350 km s ' 



detected on the jet axis, if launching occurs down to the corota- 
tion radius - 0.07 AU (Agra-Amboage et al. 201 1). Indeed, 
the predicted asymptotic speed in this model varies with anchor 
radius ro as 



^ 100 



llAU/ \0.5Mg 



1/2 



km s 



(2) 



The corresponding ejection efficiency, ^ = 0.04, also appears 
compatible with the observed ejection to accretion ratio in the 
DG Tau atomic jet, given current uncertainties (Agra-Amboage 
el al. 2011). The resulting (r, z) distribution of density and 
streamUnes is shown in Fig. 1. Thedensity contours turn from 
horizontal to vertical at polar angles < 18°, producing an ap- 
parent density collimation in good agreement with observed jet 
widths (Cabril et al. 1999; Garcia el al. 2001b; Ray el al. 2007). 
This dense axial beam is surrounded by wider streamlines from 
larger ro (white curves in Fig. 1) that recollimate on larger scale 
and reach a lower terminal speed (see eq. 2). Such a steep trans- 
verse velocity decrease from axis to edge is observed across 
several resolved Class II jets including DG Tau and CW Tau 
(Lavalley-Fouquelel al. 2000; Bacciolli el al. 2000; Coffey el al. 
2007) and seems in good agreement with the chosen MHD disk 
wind solution (Pesenti el al. 2004; Cabril 2009). 

"Slow" MHD solutions (A < 20) that become super Alfvenic 
may be obtained by including a small entropy deposition at the 
disk upper layers z - h (Casse & Ferreira 2000). The warm 
mass-loading region is denser and the ejection efficiency ^ is 
substantially enhanced, lowering the A value (see eq. 1). Such 
additional heating was assumed to be powered by dissipation of 
MHD disk turbulence and parametrized as a tiny fraction / of the 
disk accretion power, with / = 8 x 10 in the solution chosen 
here (see Casse & Ferreira 2000 for a discussion of the effect of 
/ on the wind dynamics). This new class of solutions, referred to 
as "warm" in Pesenli el al. (2004), is still "cold" in a dynamical 
sense as the initial thermal energy remains negligible compared 
to gravity, and the wind acceleration is still mostly magnetic. 
Therefore, the wind temperature beyond the slow point does not 
affect the dynamics and we may recompute it a posteriori from 
the actual heating/cooling terms, without loss of self-consistency 
(see discussions in Section 4 and Appendix A). 

Other input properties of the adopted solution are: a thermal 
(uncompressed) aspect ratio s - h/r - Cs/Vk = 0.03 (with Cj 
the sound speed and Vk the kepler speed), a turbulent resistiv- 
ity in the disk midplane Vn, - a,^VA.Qh^ with Va.o the Alfven 
velocity and - 2 (in order to ensure stationarity), and a tur- 
bulent viscosity v,, - v,„. Crossing of the slow magnetosonic 
point is obtained for a ratio of magnetic to thermal pressure in 
the disk midplane fi ^ 0.33 (leading to quasi-keplerian rotation; 
Shu el al. 2008). The computed second Blandford & Payne jet 
parameter is k ^ 0.1, and the field inclination at the disc surface 
is Oq ^ 57.3°, with similar surface values of B;, and Br. 

For a given self-similar MHD solution, all physical quanti- 
ties at a given polar angle obey specific scalings with the stellar 
mass Mi,, the disk accretion rate Mace, and the anchoring radius 
ro of the magnetic surface in the midplane, as given in eq. (9) 
of Garcia el al. (2001a). In particular, the midplane field in the 
adopted MHD solution is. 



B, ^ 0.7G 



ll AU/ 



-5/4 



Mace 



10-''Moyr- 



1/2 



0.5Mp 



^1/4 



(3) 



which follows the same scaling as the available magnetic field 
measurements in protostellar disks (Shu el al. 2007). The re- 
quired field strengths thus appear plausible. In the following, we 
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Figure 1. Overall geometry of the "slow" MHD disk wind solu- 
tion of (Casse & Ferreira 2000) used in this article. Solid white 
curves show various magnetic flow surfaces, with that anchored 
at 1 AU shown in dashed. The density for Mace = lO^^Mo/yr and 
Mi, = O.SMq, is coded in the contourplot starting at 8x 10^ cm"^ 
and increasing by factors of 2. The bottom dotted line traces the 
slow magnetosonic surface (at 1.7 disk scale heights) where our 
chemical integration starts. A color version of this figure is avail- 
able in the online edition of this journal. 

will vary M^,, Mace, and ro in order to explore their effect on the 
molecular content and thermal state of the disk wind. 

2.2. Thermo-chemical evolution 

The thermo-chemical evolution along wind streamlines was im- 
plemented by adapting the last version (Flower & Pineau des 
Forets 2003) of a code constructed to calculate the steady-state 
structure of planar molecular MHD multi-fluid shocks in inter- 
stellar clouds (Flower et al. 1985). The ion-neutral drift, FUV 
field, and dust-attenuation were calculated following the meth- 
ods developed by Garcia et al. (2001a) in the atomic disk wind 
case. Irradiation by stellar X-rays was also added, following the 
approach of Shang et al. (2002). 

Given the high densities and low ion-neutral drift speeds in 
the wind (see Section 3) the same temperature and velocity is 
adopted here for all particles. The latter is prescribed by the 
underlying single-fluid MHD wind solution. We thus integrate 
numerically the following differential equations on mass den- 
sity p, species number density «(«), and temperature T along the 



streamline, as a function of altitude z above the disk midplane: 

dp/ ^ ■^/-P/V-v 
dz V, 
dn(fl) Ra - n(a)W ■ v 
dz v- 

dT F - A - niotkeTV ■ v - ^keTR 

T' ^ V, ■ 

az f^BntotV, 

Here, V; and (V ■ v) are the bulk vertical flow speed and the (3D) 
flow divergence interpolated from the MHD wind solution, Wtot 
is the total number density of particles, T is the temperature, 
kg is the Boltzmann constant, S and R are the rates of change 
in mass and number of particles, respectively, per unit volume, 
and r and A are the heating and cooling rates per unit volume. 
The equations on n{a) apply to each species a as well as to the 
individual populations of the first 49 levels of H2 (up to an en- 
ergy of 20 000 K) which are integrated in parallel with the other 
variables. The equations on py apply to each "fluid" / (neutral, 
positive, negative) and are used mainly for internal checking pur- 
poses, as the total mass density p(z) is prescribed by the MHD 
solution. 

Cooling and heating mechanisms include 

- Radiative cooling by H2 lines excited by collisions with H, 
H2, He, and electrons (Le Bourlot et al. 1999). 

- Radiative cooling by CO, H2O, and '^CO in the Large 
Velocity Gradient approximation (Neufeld & Kaufman 
1993), and by OH and NH3 in the low-density limit (Flower 
etal. 1985). 

- Atomic cooling by fine-structure and metastable lines of C, 
N, O, S, Si, C+, N+, 0+, S+, Si+ (Floweret al. 2003) andFe+ 
(Giannini et al. 2004). 

- Inelastic scattering of electrons on H and H2 (Aggarwal et al. 
1991; Hummer 1963; Rapp & Englander-Golden 1965). 

- Energy released by collisional ionisation and dissociation 
and exo/endo-thermicity of chemical reactions (Flower et al. 
1985). 

- Energy heat/loss through thermalization with grains (Tielens 
& Hollenbach 1985). 

- Ambipolar-diffusion heating by elastic scattering between 
the neutral fluid and charged ions and grains (Garcia et al. 
2001a, see Section 2.3). 

- Ohmic heating arising from the drift between electrons and 
other fluids (Garcia et al. 2001a, see Section 2.3). 

- Photoelectric effect on grains (Bakes & Tielens 1994, eq. 42) 
irradiated by the (attenuated) FUV field of hot accretion 
spots (see Section 2.4). 

- Heating through cosmic-rays and (attenuated) coronal stellar 
X-rays (Dalgarno et al. 1999; see Section 2.5). 

The reader is referred to the corresponding references for a 
discussion of the physical context and hypotheses involved in 
modelling each of the above processes. 

2.2.1. Chemical network 

The chemical network consists of 134 species, including atoms 
and molecules (either neutral or singly ionized) as well as 
their correspondents inside grain refractory cores, and on grain 
icy mantles. A representative polycyclic aromatic hydrocarbon 
(PAH) with 54 carbon atoms is also included, with a fractional 
abundance of 10"^ per H nucleus. The total elemental abun- 
dances and their initial distribution among gas, grain cores, and 
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icy mantles are taken from Tables 1 and 4 of Flower & Pineau 
des Forets (2003). 

The charge balance of grains and PAHs is treated as in 
Flower & Pineau des Forets (2003). The number of grains per 
H nucleus and the mean grain cross section are determined 
from the abundance of depleted elements in grain cores, and 
the adopted grain size distribution (see Section 2.4.2). Most of 
the grains are charged, and assumed well-coupled to the charged 
fluid. 

We consider 1143 reactions including neutral-neutral and 
ion-neutral reactions, recombination with electrons, charge ex- 
change, cosmic ray induced desorption from grains, sputtering of 
grain icy mantles, and erosion of charged grain cores by impact 
of drifting heavy neutral species. Reaction rates between charged 
and neutral species are enhanced by ion-neutral drift, following 
the effective temperature prescription of Flower et al. (1985), 
with the drift speed v,„ computed as described in Section 2.3. 

H2 coUisional dissociation is treated level by level following 
Le Bourlot et al. (2002). Reformation on grains is computed with 
the sticking coefficient of Hollenbach & McKee (1979). Three- 
body reactions in the gas phase are not included, therefore we 
will limit ourselves to densities < lO'^ cm"-'. 

lonisation and dissociation reactions by far-ultraviolet 
(FUV) photons from accretion hot spots, and by stellar coro- 
nal X-rays, are included as described in Sections 2.4 and 2.5, 
respectively. 

2.2.2. Initial conditions and integration 

The integration of the set of equations for temperature, mass and 
number of particles along the flow is an initial value problem. 
Thus initialization of temperature and initial populations have to 
be devised. As in Garcia et al. (2001a), we start all calculations 
from the wind slow magnetosonic point (located at 1 .66h{r) for 
our adopted MHD solution), and assume that the temperature 
and species abundances there are at equilibrium with the local ra- 
diation field. These equilibrium values are obtained by perform- 
ing a "steady-state" run over 10^ yrs where thermo-chemical 
equations are solved with the density held fixed. We then in- 
tegrate the thermo-chemical equations along the flow using the 
DVODE integrator (Brown et al. 1989), until the recollimation 
point where the streamline reaches its maximum radial extension 
(radius of 40ro at z = 900ro for our chosen wind solution). We 
checked that the final temperature and H2 abundance along the 
streamline do not depend sensitively on the initial equilibrium 
conditions, as long as the gas is fully molecular initially. This oc- 
curs because the ionisation (which controls the ion-neutral drift 
heating) adjusts rapidly at the dense wind base to the local radi- 
ation field, even though it becomes frozen-in at large distances. 



2.3. Ambipolar diffusion coefficients 

The ambipolar diffusion and Ohmic heating terms require a de- 
tailed calculation of elastic momentum exchange rates and drift 
speeds. The ion-neutral drift speed is calculated from the MHD 
solution using the analytical formula derived in Appendix A of 
Garcia et al. (2001a), valid for ||v,„|| « ||v||: 



1 



JxB 



1 + p \c{\ + X)fj.i„niVi, 



+ P- 



qeth 



p = 



(7) 
(8) 



where J is the electric current density, B is the magnetic field, c 
is the speed of light, q,, is the electron charge, rig and n, are the 
number densities of electrons and ions, X = Pi/pn, and ju,„«,v,„ 
and fifniigVen are the total momentum transfer rate coefficients for 
ion-neutral and electron-neutral collisions, respectively. 

The total ion-neutral momentum transfer rate coefficient is 
obtained by summing over the main neutrals (H, H2, He) and 
over all ions, charged PAHs and charged grains: 

^ ^ jJahn{a)n(b){o-v)„i,, (9) 

a - ions, b = 
PAH^,^^ H,H2,He 



PinniVin = 



where n{a) is the number density of species a, pai, is the reduced 
mass of particles a and b, and {o-v}gi, is their elastic collision 
rate coefficient; the latter is evaluated via the recent analytical 
fits to quantum-mechanical calculations provided in Table 2 of 
Pinto & Gain (2008a), when available, given as a function of the 
effective relative speed v: 



V - 



+ \n 



with Tr 



T^Pah 



(10) 



(11) 



For the H2-H^ pair we use the updated fit in Pinto & Galli 
(2008b). For the rest of the pairs {crv}„i, is taken as the maximum 
of the polarizability and hard sphere rate coefficients (Garcia 
etal. 2001a): 



max \ l.Alnqe — , ctsV 

Pab 



(12) 



where cc/, is the polarizability of H, He, or H2 and the hard sphere 
cross section is taken as 



10-"* cm^ , PAH^ 
^{rI) , grains (g^) 
10"''' cm- , else 



(13) 



Here ^(Rj) denotes the average grain cross section over our 
size distribution. The total momentum transfer rate coefficient 
for collisions between electrons and neutrals is given by: 



PenneVen = rUeHe ^ «(a) {(Tv) ^ 
a=H,HT,He 



(14) 



with the {(Tv)ga also evaluated from the recent analytical fits of 
Pinto & Galli (2008a). The ion-electron momentum transfer rate 
entering Ohmic heating is calculated with the classical formula 
of Schunk (1975) (see eq. A.13 of Garcia et al. 2001a). 

Fig. 2 plots the momentum transfer rate coefficients (nor- 
malized by n^) for collisions between the neutral fluid and the 
different charged fluids as a function of vertical distance above 
the disk midplane for a representative streamline. It may be seen 
that collisions with positive ions dominate the coupling, with an 
extra contribution from charged PAHs near the flow base. 

The ambipolar drag heating and Ohmic heating are given by 
(Garcia etal. 2001a): 



^ drag — 1 



2 1 



(15) 



with T] = (pieiiVie + pe„n„v„t.) I {q,,n^f . Note that fdrag is pro- 
portional to the square of the Lorentz force per particle, and in- 
versely proportional to the ionisation fraction in the wind. 
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Figure 2. Momentum transfer rate coefficients (normalized by 
in collisions of the neutral fluid with the different charged 
fluids as a function of altitude along a flow line launched at 1 AU, 
for Mace = 10-^Mo yr-\ = Q.SM^ (our Class 1 model). A 
color version of this figure is available in the online edition of 
this journal. 



2.4. FUV field 

The FUV radiation will control the ion abundance (and therefore 
the ion-neutral drift) in the upper wind regions. It will also pho- 
todissociate molecules, and provide extra heating by collisions 
with warm dust and photoelectric effect on grains. Our treat- 
ment essentially follows the approach of Garcia et al. (2001a) 
with several updates, including an accretion spot geometry and 
a photochemical network. 



2.4.1. Unattenuated radiation field 

Instead of the optically thick equatorial boundary layer model 
adopted in Safier (1993) and Garcia et al. (2001a), we adopt here 
the paradigm of magnetospheric accretion, recently favored in 
T Tauri stars, where the inner disk is truncated at several stel- 
lar radii and accreting material flows along field lines in quasi 
free-fall, giving rise to shocked "hot spots" on the stellar surface 
(Bouvier et al. 2007). The accretion geometry in younger pro- 
tostars is not yet constrained observationally, but their similar 
X-ray flare properties to T Tauri stars indicate an akeady well- 
developed stellar magnetosphere (Imanishi et al. 2003), hence 
we will adopt the same paradigm for simplicity. 

The FUV radiation field thus comes from the central low- 
mass star of effective temperature = 4 000 K, and from hot 
accretion spots of fixed temperature Ths - 10000 K, chosen to 
match the mean observed colour temperature of the FUV excess 
in T Tauri stars (Johns-KniU et al. 2000). Both radiation sources 
are treated as black bodies. 

As in Garcia et al. (2001a), we neglect the scattering contri- 
bution to the FUV field^. The (unattenuated) direct stellar flux at 
distance R from the star (with dimensions erg cm"^ s~' Hz~') is 
calculated through the relation 



(16) 



- Radiative transfer calculations in 1 + lD show that it starts being 
important only very near the disk plane (Nomura & Millar 2005). 
Ionization and dissociation in this region will be dominated by hard 
stellar X-rays in our model. 



where By(T) denotes the specific intensity at frequency v of a 
black body of temperature T, Ri, is the stellar radius, and 5hs is 
the fraction of the total stellar surface 4-jtRI. covered by hot spots 
(assumed uniformly distributed). The main difference with the 
boundary-layer model of Garcia et al. (200 1 a) is thus a fixed Ths 
independent of Mace, and a more isotropic radiation flux (due to 
the lack of projection effects or disk occultation). 

In each model, 6hs is determined by requiring that 
the hot spots radiate half of the accretion luminosity, 
i.e. 4jTRl6h,o-BT^^ = GM^^^M-^/ilR-^). The value of one half is 
meant to be illustrative only: the actual fraction of accretion lu- 
minosity radiated in the accretion shock could exceed 80% for a 
large disk truncation radius > 57?*, or decrease if a sizeable frac- 
tion of accretion energy is tapped to drive a stellar wind. None 
of these effects being well quantified, especially in Class and 
Class I sources, we adopt 50% here for illustration and consis- 
tency with the boundary-layer model of Garcia el al. (2001a). 
We set Ri, = 3Rq, the radius of a young accreting 0.5Mo star 
near the "birthline" with T*, = 4 000 K according to the models 
of Stahler (1988), noting that it remains a good approximation 
for a very young Class (/?* - 2Rq for a O.IM© protostar ac- 
creting at 5 X lO"''M0yr"'; cf. Stahler 1988). With the above as- 
sumptions, the modelled FUV hot spot continuum in the domain 
1 400 - 2 000 A is a factor 1 .5 - 2.5 stronger than the accretion 
shock models of GuUbring et al. (2000) for DR Tau and DG Tau, 
using the same accretion rates (Mace = 3 - 5 x 10"^ Mq yr"'). 
On the other hand, if the FUV spectrum below 1 400 A were as 
flat as in the lower accretion stars BP Tau and TW Hya (Bergin 
et al. 2003), our blackbody model would underestimate the flux 
at 1 000 A by a factor 4-2. This remains sufficiently accurate 
for the present exploratory study. 

Since we do not solve for radiative transfer, we neglect the 
extra FUV flux in the stellar Lyo' line. Lya pumping of highly 
excited levels of H2 is occuring close to the star (Herczeg et al. 
2002, 2006; Nomura & Millar 2005) but will be less important 
on the larger jet scales of interest here, especially since we focus 
on dense molecular jets where H2 is well shielded. 



2.4.2. Attenuation by dust 

The stellar FUV field is attenuated mainly by dust in the disk 
wind along the line of sight to the star The sublimation radius 
R^ub is calculated via eq. (B.5) of Garcia et al. (2001a), consider- 
ing a unique sublimation temperature Tjub = 1 500 K. We used 
the dust optical properties tabulated by Draine & Lee (1984); 
Draine & Malhotra (1993); Laor & Draine (1993), assuming the 
standard "MRN" mixture of astronomical silicates and graphite 
and the grain radius distribution n{d) oc o^^^ proposed by Mathis 
et al. (1977), with grain radii in the interval [0.005,0.25] jum. The 
sublimation surface for the Class O/I/II models in Section 3 is 
plotted in Fig. 3 on top of inner disk wind streamlines. 

The dust optical depth to the star at point of spherical coor- 
dinates {R, 6) is TyiR, 9) = o-d(v)N'^{R, 6) where cr^(v) is the dust 
extinction cross section per H nucleus (integrated over the grain 
size distribution), and N^{R, ff) is the column density of hydro- 
gen nuclei on dusty wind streamlines along the line of sight to 
the star In our self-similar wind model, the wind density varies 
as R^^-^f{0), and N^{R, 6) may then be integrated analytically as; 



N^^{R, 0) = 2«h(^, e)R 



R 



1 



(17) 



Here, riu is the number density of H nuclei in all forms, 
riii ^ n(H) + 2n{H2) + n(H^), and Rminiff) is the spherical radius 
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Figure 3. The dust sublimation surface in our Class O/I models 
(outer red circle) and our Class II model (inner green circle), for 
the grain size distribution and mixture of silicates and graphite 
proposed by Mathis et al. (1977). The wind slow magnetosonic 
surface (denoted SMS) is superimposed as well as three flow 
lines (orange dotted lines), including the innermost one launched 
from 0.07 AU. The stellar parameters are = 4000 K and 
Ri, = 3^o, while the hot spot temperature and coverage fraction 
are Ths = 10000 K, and 6h, = 0.032 (Class O/I) and 0.0032 
(Class II). A color version of this figure is available in the on- 
line edition of this journal. 



inside which there is no dust, i.e. either the sublimation radius, 
Rsah, or the radius of the innermost flow line at angle 0, Rin{0), 
whichever is larger (see Fig. 3). The innermost wind streamline 
is assumed launched from ri„ - 0.07 AU, a typical corotation 
radius in young low-mass stars. 

The total dust geometrical cross section per H nucleus, cal- 
culated from our grain size distribution and the abundance of 
depleted elements in grain cores, is n^R^'^rid/nu - 1.36 x 
10"^' cm^, where (ifj^ is the mean square grain radius; the ratio 
of Ay to A^H is then cr^(5500 A) = 2.72 x lO"^^ ^.^2 

2.4.3. Dust temperature 

Collisions with dust grains heated by the UV field can be an im- 
portant gas heating/cooling term at the wind base (though unim- 
portant further out). We assume that above the wind slow point 
(z > 1 Jh), dust is optically thin to its own radiation. The dust 
temperature is then calculated at every step from the grain radi- 
ation equilibrium against incident (wind-attenuated) stellar pho- 
tons: 



AttctbT', {r^) ■ = j trrf(y)F,(^)e-<'(''><(«'«'dy 

(18) 

where (Tb is the Stefan-Boltzmann constant, (r^^ is the mean 

square grain radius, {QemiTj)) is the Planck-averaged grain 
emission efficiency (weighted by n{a)a^), and cr^(v) is the dust 
absorption cross section per H nucleus defined earlier. Collisions 
with the gas are neglected in the grain thermal balance, as they 
compete with radiation only at densities > lO'"' cm"-' (Glassgold 
et al. 2004). We verified that when the FUV excess and the wind 
attenuation are negligible, our values agree with detailed cal- 



culations at the top of the disk atmosphere by D'Alessio et al. 
(1999) for MRN dust and a 4000 K star (taking into account 
differences in i?*). 

2.4.4. Photochemical network 

A network of photoionisation and photodissociation reactions is 
implemented with rates in the form kp - xj^^^'^''\ where Ay 
is the visual extinction calculated above (Section 2.4.2), p, y are 
constants depending on the reaction, taken from van Dishoeck 
(1988) and Roberge et al. (1991), and x is the ratio at 1 000 A of 
the unattenuated FUV flux to the mean interstellar radiation field 
of Draine (1978), integrated over 4;r(2 X 10"^ ergs"' cm"^A"'). 
Such a "monochromatic" flux normalization at 1 000 A was cho- 
sen rather than the average ratio to the interstellar field over 
the interval 910 - 2 066 A, denoted Go in the literature (see 
e.g. Tielens & Hollenbach 1985), as photoreactions contributing 
to the main ionisation and dissociation reactions in our models 
(C, S, CH+, H2, CO) occur in the narrow range 910- 1 200 A. 

Special cases of photodissociation reactions are the dissoci- 
ation of H2 and CO, which occur through line absorption. We 
thus need to consider not only the shielding caused by the dust 
particles, but also the additional self and cross shielding due to 
those species themselves. They were implemented in the model 
in an approximate way as described below. 

2.4.5. H2 photodissociation. 

We adopted an H2 photodissociation rate coefficient per 
molecule of the form proposed by Draine & Bertoldi (1996): 



^P,H2 ^XPoi^i) X/h,,j(Av.)/h,,h,(A^H,,^) 



(19) 



where /:'o(H2) - 4.2 x 10 " s ' is the unshielded rate in the 
Draine field, fn^jiAy) - t^l^''^" is the shielding factor by dust, 
with Pd - crrf(1000 A)/cr,/(5500 A) = 6.3 for our dust mixture, 
and /h2,H2(/Vh2> ^) is the H2 self-shielding factor given by eq. (37) 
in Draine & Bertoldi (1996), with the line width parameter b set 
here to the local sound speed Cj. The shielding column of H2 
molecules on the line of sight to the star, A^h, , depends on the H2 
abundance on all inner wind streamlines. In order to allow a local 
calculation in the present exploratory study, we assume a mean 
H2 abundance on the line of sight to the star equal to half the 
local H2 abundance at the current point of the streamline. This 
"local" assumption tends to overestimate self-shielding, but will 
not affect our results when photodissociation is not the dominant 
destruction mechanism. A conservative check of this hypothesis 
is made a posteriori (see Section 3). 

2.4.6. CO photodissociation. 

For CO, one needs to consider not only the self-shielding and 
the shielding by dust, but also the shielding by H2 (through line 
overlap with CO). We have fitted the data of Table 1 1 from Lee 
et al. (1996) as a function of Ay and of the shielding columns 
A^co and A^H2. ™d obtained the following analytical formulae for 
the shielding factors by dust, CO, and H2 respectively: 



/co,d(Av) = max{0.5e"^-5'^",0.1e"3'*''}, 
/co,co(/Vco) = min{l,(^^^) }, 
fmm.) = O.Se- 



dOi5 cm-2, 

-A'h,/1.7x10-' cm"- 



(20) 
(21) 
(22) 
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Figure 4. Adopted analytical fits (red curves, from eq. 22) to the CO photodissociation shielding factors of Lee et al. (1996) (con- 
nected green symbols) caused by CO line overlap (left), H2 line overlap (middle) and dust extinction (right). 



These analytical formulas are valid to within a factor of 2 up to 
A^co - 2 X 10''' cm"-', and over the whole range in A^h, and Ay, 
as shown in Figure 4. 

The photodissociation rate coefficient per molecule is then 

^p.co = XPo(CO) X /co,rf(A v)/co,co(A^co)/co,H2 (A^h, ), (23) 

where /?o(CO) = 2 x 10""^s"' is the unshielded rate for an 
isotropic Draine field (twice the value from Lee el al. (1996), 
who assume a one-sided cloud illumination). As done for H2, 
the shielding column of CO molecules on the line of sight to the 
star, A^co, is calculated assuming a mean CO abundance on inner 
disk streamlines of half the local value of the CO abundance at 
the current point. 

2.5. X-rays 

Hard coronal X-rays will be the dominant ionization process at 
the base of the dusty disk wind, where the stellar FUV flux is 
strongly extinguished by dense inner streamlines. Energetic sec- 
ondary electrons generated by X-ray ionisation will also dissoci- 
ate and heat the gas. Our adopted X-ray treatment and chemical 
network are described in the following. 




10"^ 10'^ 10"'' 10° 10^ 10^ 10^ 



Figures. The function Jp(Tx) describing the attenuation of the 
energy deposition rate by a thermal X-ray spectrum versus tx, 
the X-ray optical depth at the characteristic energy, kTx- Also 
plotted is the function TxJpirx) describing the effect of chang- 
ing kTx (i-e- Tx) for a fixed attenuating gas column density (see 
eq. 28). 



where 



2.5. L X-ray spectrum and attenuation 

The X-ray flux is modelled as a thermal spectrum with charac- 
teristic energy kTx, 



Ho(R) 

Ipi.^O,Tx) 



70-r,e(kTx), 



r 



(26) 
(27) 



Fx(E) = dLx{E)/dE = (Lx/kTx) x exp(-E/kTx) 

The attenuated rate of X-ray energy deposition per H nucleus in 
the wind writes 

Hx{R, B) = ^^Pe(£) exp-p=^^>^"<«'«' d£, (24) 

where Eq is the low energy cutoff of the incoming X-ray spec- 
trum, crpg(£') is the photoelectric cross section per H nucleus, 
and Nh{R, d) is the total column of H nuclei to the star through 
inner wind streamlines (given by eq. 17 with RmmiO) - Rin{0)). 
We calculate Hx analytically following Glassgold et al. (1997) 
and Shang et al. (2002), who adopted a low-energy cutoff 
£0 = 0.1 keV and a power-law approximation to the cross sec- 
tion, crpe(£) = 2.27 X 10-22 cm-2 (KeV/£)'' with p = 2.485. Hx 
then writes 

Hx(R,6) = HQ(R)Ip{^o,Tx) (25) 



^0 = Eo/kTx, Tx = ap^(kTx)Nu{R,0) is the X-ray optical 
depth at kTx, and Jp{Tx) is an analytical fit given in eq. (C3) 
of Shang et al. 2002, and plotted in Fig. 5. The fit is valid 
for Tx> tq - 1/p^q{p + ^0) = 10"^, a condition verified in our 
models. 

The effect of changing kTx for a fixed Lx and fixed posi- 
tion (R, 6) (and therefore a fixed attenuating column A^h) may be 
visualised by rewriting Hx as 

HxiR, 6) = Ho{R)Jp{Tx) = TxJpiTx) x Lx/i^nR^Na). (28) 

The function TxJp{Tx), also plotted in Fig. 5, is almost flat over 
a broad range of tx (0.001 to 100). The regions where X-ray 
ionisation dominates in our models have values of tx - 0.3 - 30 
within this range for kTx - 4 keV. Furthermore, the ionization 
fraction varies only with the square root of Hx. Varying kTx in 
the range 2.5 - 15 keV (i.e. tx in the range 0.01 - 100) would 
thus not strongly modify our results. 
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2.5.2. X-ray ionization rates 

X-ray ionization rates are calculated with the same simplify- 
ing assumptions as in Glassgold et al. (1997) and Shang et al. 
(2002): The mean energy to form an A* ion is approximated by 
its constant high-energy limit W(A'^), valid for primary photo- 
electrons with E > Eq - Q.l keV. Neglecting the small differ- 
ence in energy between the absorbed X-ray photon and the pri- 
mary photoelectron, and recalling that the absorbed energy Hx 
is per H nucleus, the X-ray ionization rate of species A per atom 
is then given by 

HyHu 1 

The limiting mean energies W(j4^) for X-ray ionization of 
H2, H, and He in neutral gas mixtures with 10% of Helium have 
been calculated by Dalgarno et al. (1999) who fitted them"* as 
linear functions of «(H)/n(H2) (middle terms of eq. 30-32 be- 
low). We further simplify their expressions by taking 0.53 - 0.5, 
1.89 - 2, and recalling that the number density of H nuclei 
hh - «(H) + 2«(H2) at low ionization, and «(He) = 0.1«//, which 
yields the right-hand side terms: 



W(Y{^_) = 41.9eV 1 -H 0.53 

WiW^) = 39.5eV 1 + 1.89 
H^(He+) = 470eV ^ 47.0eV 



«(H) 



n(H2) 

«(H2) 



«(H) 
n(He)' 



^ 20eV 
^40eV 



«H 
«(H2)' 

«(H)' 



(30) 

(31) 
(32) 



Inserting the right-hand side expressions for W(A^) in 
eq. (29), the X-ray ionization rates per specie A become inde- 
pendent of «H and n{A) and may be written solely in terms of the 
parameter ^(Hj) that enters our chemical network: 



an^) - 0.5^(H+), 
^(He+) - 0.4f(H+). 



(33) 
(34) 
(35) 



The total ionization rate per H nucleus, often denoted in the 
literature, is related to ^(H:^^) by 



- mi) 



0.5«(H) + «(H2) + 0.4«(He) 



0.54f(H^), 



(36) 
(37) 



regardless of the ratio H/H2. The formation of H^ through disso- 
ciative ionization of H2 is also included with a rate of ^{H^)/22 
per molecule (Dalgarno et al. 1999). A constant cosmic ray ion- 
ization rate fcR(H^) = 5x10"'^ s"' is added to the rate produced 
by stellar X-rays, but it plays a negligible role in our models. 



2.5.3. X-ray induced dissociation 

The rate of H2 dissociation by X-ray induced electrons is taken 
as 0.5^(H+), following the results of Dalgarno et al. (1999) for 
primary electron energies > 0.1 keV. The energetic electrons 
also collisionally excite H2 to higher electronic levels, which 
then radiatively decay by emitting a flux of ultraviolet "sec- 
ondary photons" able to dissociate other molecular species. We 
adopt the dissociation factors pM calculated by Gredel et al. 



' The variation of IV(A^) with ionization fraction is neglected here, 
as it amounts to less than 3% for the ionization fractions < 10"' en- 
countered in our models (Dalgarno et al. 1999) 



(1989) for various species, leading to a rate per unit volume of 
^(H+)[«(H2)/«H]n(M)/7M/(l - w) (Flower et al. 2007), where 
we adopt a grain albedo ai = 0.5. For CO, we adopt the high- 
abundance, high-temperature limit of pco - 28 for co - 0.5 
(Gredel et al. 1987). Photodetachment of electrons from PAHs 
by secondary photons is also included (Flower & Pineau des 
Forets 2003). 



2.5.4. X-ray lieating 

The heating efliciency of the X-rays and cosmic rays is treated 
collectively. We use the results of Dalgarno et al. (1999), who 
studied the heating produced by high energy electrons interact- 
ing with a partly ionized gas mixture of H, H2 and 10% of He. 
The heating efficiency ///, is calculated for ionization fractions 
< 0. 1 and includes the fraction of initial primary electron energy 
lost by elastic scattering with neutrals and electrons and by rota- 
tional excitation of H2 (rapidly thermalized at our densities). The 
total heating rate per unit volume by X-rays and cosmic rays is 
given by: 

lOrT/H, + rjii 

Tx+cr(z) = TjhHxnu = — - — f(H,)eH,«H (38) 
lOr +1 " ' 

where e^^ - 20 eV follows from the definition of ^(Hj) in 
eq. (33), r — «(H2)/«(H), and //h,, '7h are respectively the heating 
efficiencies for the H2-He and the H-He ionized mixtures, fitted 
as a function of the electron fraction .^e = «e/nH in the form 



1 + 



?7o ■ 



1 



1 + cxl 



(39) 



with 770 = 0.055, c = 2.170, a = 0.366 for H2-He, and 
770 = 0.117, c = 7.950, a = 0.678 for H-He (Dalgarno et al. 
1999, Table 7). 



3. Results 

To investigate how the chemical content of protostellar disk 
winds evolves in time with the decline of accretion rate, we con- 
sidered three sets of parameters representative of the main evo- 
lutionary stages of young low-mass stars with bright jets: 

- a very young Class protostar (Mace = 5 x lO^^'M© yr"'. 
Mi, - 0.1 M0), where accretion is high and the star has not 
yet reached its final mass; the chosen parameter values de- 
scribe e.g. the heavily embedded sub-millimeter exciting 
source of the HH 212 jet (Lee et al. 2006). 

- a Class I source (Mace = 10-'' Mo yr-i, M* = 0.5Mo), where 
the star has accumulated most of its final mass, but resid- 
ual infall and accretion still proceed at a fast pace; the cho- 
sen parameters describe e.g. the infrared exciting source of 
HH26 (Antoniucci et al. 2008). 

- an active Class II star (Mace = IO^^Mq yr"', M* = O.5M0), 
representative of optically visible T Tauri stars with bright 
optical microjets, e.g. DG Tau or RW Aur (GuUbring et al. 
2000). 



With our choice of parameters, and the self-similar scalings of 
the MHD solution, the wind density «h oc MaccM^'^^r^ drops 
by a factor 10 from one Class to the next. Our adopted stellar 
radius of 3/?© gives the same accretion luminosity of 5Lo in the 
Class and Class I models, and 10 times smaller in the Class II, 
leading to a hot spot coverage fraction (5hs of 3.2% and 0.32% 
respectively (cf. Section 2.4.1). 
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Figure 6. Various variables along a flow line anchored at 1 AU, for the dense MHD disk wind solution and three typical classes of 
stars: class (M* = OAMq, M-^^c - 5 x IO'^Mq yr"'; dotted blue curves), class 1 (M* = O.SMq, Mace = IO^^Mq yr"'; solid red 
curves) and active class II (M* - Q.5Mq, Mace = IO^^Mq yr"^; dashed green curves), (a) The gas temperature T; (b) the number 
density of hydrogen nuclei «h; (c) the visual extinction Ay to the star, in magnitudes; (d) the poloidal components of the bulk flow 
speed V (smooth top curves) and calculated drift speed v,„ (lower 3 curves); Note that the Class II and Class I have the same M*, 
hence the same bulk flow speed for ro = 1 AU; (e) the X-ray H2 ionization rate per molecule f (Hj ); (f) the effective radiation field 
in Draine units ;(fexp"^'^'' (attenuation factor appropriate for carbon photoionisation); (g) the fractional abundance of ions; (h) the 
fractional abundance of H2 (upper limit in the Class II model, see Section 3. 1.3); (i) the fractional abundance of CO (upper limits in 
the Class I/II models, see Section 3. 1 .4). A color version of this figure is available in the on-line edition of this journal. 



We adopt a thermal X-ray spectrum of characteristic energy 
kTx = 4 keV and luminosity Lx = 10^'^ erg/s, typical of the 
"hard" time-variable component in solar-mass young stars and 
jet-driving protostars (Imanishi et al. 2003; Gudel et al. 2007). 
Lx is a time-averaged value of 30% of the median flare luminos- 
ity in solar-mass YSOs, based on a typical flare interval of 4-6 
days (Wolk et al. 2005). The exact value of kTx has no major 
effect on our results (see Section 2.5.1). 



Note that the flow crossing timescale out to the recollimation 
point is very short, only 50 yrs for a streamline launched at 1 AU 
in the Class I model. The chemical composition of the disk wind 
will thus deviate substantially from published "static" disk at- 
mosphere models such as, e.g. Glassgold el al. (2004); NomLira 
& Millar (2005). The presence of adiabatic cooling and drag 
heating also introduce major differences in the thermal structure. 
These effects are discussed in Section 4. 



We first analyse in detail the results obtained for the 3 evo- 
lutionary classes, for a representative anchor'* radius ro = 1 AU, 
and for a streamline launched just outside sublimation, ro - Rsub- 
The eff'ect of increasing ro is discussed later on, in Section 3.2. 

3. 1 . Results for ro - I AU and ro = ^sub for the 3 classes 

Figure 6 shows the main physical parameters integrated as a 
function of z along the 1 AU streamline for the Class O/I/II mod- 
els: temperature, ion-neutral drift speed, ionisation fraction, H2 
and CO abundance. Also plotted are the (prescribed) flow den- 
sity and bulk speed, and the calculated Ay to the central star, 
attenuated X-ray ionisation rate of H2, and "effective" FUV field 
at 1 000 A in Draine units (e"-'''*'' being the attenuation factor for 



The anchor radius rg is the radius of the magnetic field surface in 
the disk midplane. The launch radius at the slow magnetosonic point 
is 6.5% larger — allowing centrifugal acceleration of the matter loaded 
onto the field line. 
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Figure?. Same as Figure 6 for streamlines anchored just beyond the sublimation radius: ro = ^sub = 0.2 AU in the Class and 
Class I, and 0.11 AU in the Class II. These calculations assume no self-shielding of H2 or CO (i.e. no molecules surviving on 
streamlines launched inside /?sub)- The resulting abundances set a minimum photodissociation timescale on outer streamlines, to be 
compared with the flow timescale (see Figure 10). A color version of this figure is available in the on-line edition of this journal. 



Carbon photoionisation in our chemical network). In Figure 7, 
we show for comparison the results obtained for streamlines 
with an anchor radius just beyond the sublimation radius, com- 
puted assuming no H2 or CO self-shielding (i.e. assuming no 
molecules on streamlines launched inside 7?sub)- 



3.1 .1 . Radiation field and ionization structure 

The ionization fraction is a key parameter governing the tem- 
perature profile, through the drag heating term. It is controlled 
by attenuation of the radiation field through inner wind stream- 
lines. The behavior is qualitatively similar along the streamlines 
launched from 1 AU and from TJsub- 

The Ky towards the star is highest at the wind base, where 
the line of sight crosses the densest wind layers (see Fig. 1) 
and drops steadily as material climbs along the streamline 
(Figure 6-7c). This causes a steep rise in the effective at- 
tenuated radiation field, until the l//?~ dilution factor takes 
over (Fig. 6-7e,f). Combined with the wind density fall-off 
(Fig. 6-7b), this generates a global increase in ionisation fraction 
along the streamlines out to 30 - 100 AU, until recombination 
sets in (Fig. 6-7g). In our models, ionization is dominated by 
stellar X-rays for Ay > 3 mag, and by FUV photons further out. 



As the wind density drops by a factor 10 from Class 
to Class I to Class II (Fig. 6-7b), the ionisation fraction, 
X{\^) = n(ions)/nH, increases by a factor 3 to 10. Indeed, the 
smaller attenuating column through the inner wind (Fig. 6-7c) 
leads to higher X-ray and FUV ionization rates. At the same 
time, the lower wind density reduces recombination. Both ef- 
fects work in the same direction. 

Fig. 8 shows the main ionization contributors along the 1 AU 
streamline for the Class 0, Class I, and Class II models. From 
these results, it can be seen that the main contributors depend on 
the overall ionization degree, X{\^): 

when X(i^) < 10 The main contributors to the negative 
charge are the PAHs. The main positive carriers are mole- 
cular ions, with the most abundant being HSO^ and HsO^; 
and (a direct product of X-ray ionization) also contribut- 
ing in the Class II. 

when X(i^) > 10"^: The main contributors to the negative 
charge are the electrons. The main positive carriers are the 
atomic ions S^, C^, and H^, which recombine more slowly 
with electrons than molecular ions do. 

An interesting result of our calculations is the relatively high 
abundances of CH+, SH+, and H+ of 10"^ - 10 reached along 
the Class I and Class II streamlines (Fig. 8). CH^ and SH^ are 
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z (AU) z (AU) 

Figure 8. The main ionization contributors along a disk wind streamline launched from 1 AU, for our Class 0, Class 1, and Class II 
models (top to bottom). The left panels plot the fractional abundances relative to H nuclei of molecular ions; the right-hand panels 
plot abundances of atomic ions and negative species (electrons and PAH ). The total ionisation fraction, X{i^) = n(ions)/«H, is also 
plotted as a thick black curve. A color version of this figure is available in the on-line edition of this journal. 



formed by endothermic reaction of C^ and with H2, when the 
wind temperature exceeds about 3 000 K. forms by X-ray 
ionisation and charge exchange, but also by photodissociation 
of CH^, with carbon playing the role of a catalyst through the 
reaction chain: 

C + hv ^ C+, (40) 
C-" H- H2 ^ CH+ H- H - 4640 K, (41) 
CH^ +hv C-HH+. (42) 

Through this repeated formation cycle, H2 is partly converted 
into H^, helping to offset destruction by charge exchange. 
The flow crossing timescale is also too short for significant 
radiative recombination. This example clearly illustrates the out- 
of-equilibrium and hybrid nature of the chemistry in MHD disk 
winds, combining elements from both C-shocks (strong heating 
by ion-neutral drift) and photo-dissociation regions (abundant 
C+ and S+). 



3.1.2. Drift speed, heating/cooling processes, and 
temperature behavior 

Thanks to X-ray and UV ionization, we find that the drift speed 
remains relatively low along the streamlines, at 10% of the bulk 
poloidal speed or less, as seen in Fig. 6-7d. This validates a 
posteriori the single fluid approximation made in computing the 
MHD dynamical solution. It also keeps the drag heating and gas 
temperature to a moderate value. Note that the ratio of drift speed 
to poloidal speed scales as Vi„/Vp oc M*M^j,|,X(i^)"' < crv >7^', 
so that the single fluid approximation remains fulfilled in the 
Class jet (small M*, large Mace) despite its low ionisation. 

The main heating and cooling terms for ro - 1 AU are plot- 
ted in Fig. 9. Collisions with warm dust can be important near the 
wind base when densities exceed a few 10"^ cm"-', but ambipolar 
"drag" heating by ion-neutral collisions, Fjiag, quickly prevails 
as density drops. The predominance of ambipolar diffusion over 
other sources of heating (shocks excepted) is a widespread char- 
acteristic of steady MHD-driven protostellar winds from low- 
mass sources, and results from the strong accelerating J xB force 
combined with a low wind ionisation (see Ruden et al. 1990; 
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Figure 9. The main cooling (left panels) and heating (right panels) processes along a disk wind streamline launched from 1 AU, 
for our Class 0, Class I, and Class II models (top to bottom). The heat transfer from grain to gas, Fg,,,, appears in the cooling 
curves when negative, and in the heating curves when positive. Other plotted cooling/heating terms are: Aadia^ adiabatic cooling; 
Ah2,h,o,co: line cooling from H2, H2O, CO; Fdrae: ambipolar diffusion heating; Fohm: ohmic heating; Fg photo^ photoelectric effect on 
grains; Fx.ci- X-ray and cosmic ray heating. OH line cooling never dominates and is not plotted for clarity. A color version of this 
figure is available in the on-line edition of this journal. 



Safier 1993; Garcia et al. 2001a). Cooling is dominated by adia- 
batic expansion out to z ^ 1 AU, and by radiative cooling further 
out, due to the high molecular abundances. The main coolant is 
H2 in the Class I and Class II models, and H2O in the Class 
model. 

The resulting behavior of temperature along the 1 AU 
streamlines is qualitatively similar regardless of the accretion 
rate (Fig. 6a): After initial cooling due to adiabatic expansion, 
drag heating takes over and heats up the gas to 1 000 - 5 000 K. 
This leads to a sharp increase in molecular cooling (denoted Aiad 
thereafter) which (assisted by collisional dissociation cooling in 
the Class I/II models) eventually overcomes Fdiag, leading to 
a temperature turnover followed by a slow decline where both 
terms remain in close balance. 

The balance between Fdiag and A,ad on the 1 AU wind 
streamline is maintained thanks to the powerful thermostatic 
behavior of molecular cooling, which is highly sensitive to 
temperature. Any slight excess/deficit of drag heating leads 
to an increase/decrease in radiative cooling until the two 



terms match again. In the outer wind, radiative cooling ap- 
proaches the low-density limit (i.e. Arad - «^/rad(r)), there- 
fore the asymptotic jet temperature profile is determined by 
the condition /,-ad(7") - Tdi-ag/ny. The latter term slowly de- 
clines with distance in our models, leading to the slowly de- 
clining temperature. Furthermore, the self-similar scalings for 
/, B, «H in eq. (9) of Garcia et al. (2001a) yield a scaling for 
Tdrag/n^ ^ ■'^*-'^acc'"o'^(i^) explaining why we find much 
cooler jets in the Class case. Final temperatures are only 
~ 500 - 700 K, compared to 2 000 K in the Class I, and 3 000 K 
in the Class II models. 

The streamlines from ro = Rsuh in the Class I and II 
cases show a slightly different asymptotic temperature behavior, 
with a flatter isothermal "plateau" at large distances (Fig. 7a). 
Molecules are heavily photodissociated on these two stream- 
lines, so that ambipolar diffusion heating is balanced by adia- 
batic cooling alone. The flat temperature "plateau" resulting 
from this balance between Fdrag and Aadia is a well known typ- 
ical property of atomic self-similar MHD disk winds (Safier 
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z(AU) 

Figure 10. The net timescales (after subtraction of any re- 
verse reaction rate) of the principle reactions that form or 
destruct H2 are plotted versus altitude along a streamline 
launched at 1 AU for our Class I model. The flow cross- 
ing timescale ff is also plotted at the bottom for comparison. 
Thin curves denote reactions whose net effect leads to destruc- 
tion, and thick curves those that lead to net creation. The rea- 
ctions are H -1- H = H2: catalytic formation of H2 on grains; 
H2 -F hv = H H- H: photodissociation (with the minimum self- 
shielded rate shown on a separate curve); H2 + x = H2 + e : 
X-ray ionisation; O + 2H2 = H2O + 2H: combined oxydation 
reactions with O and OH; H2 -H H and H2 + H2: collisional disso- 
ciation. The latter becomes very efficient near z - 2 AU, where 
the wind temperature reaches 4 000 K. A color version of this 
figure is available in the on-line edition of this journal. 



1993; Garcia et al. 2001a). The lower atomic plateau temper- 
ature ^ 4000 K found here, compared to ^ 10"* K in Garcia 
et al. (200 1 a), results from our different adopted MHD wind so- 
lution, which is typically 2 times slower and 5 times denser (see 
Section 2.1). 

It is also noteworthy that, for both values of ro, the initial 
temperature rise is much slower in the Class jet; as shown in 
Garcia el al. (200 1 : ) (equations 24 and 26), the slope of the ini- 
tial temperature rise, where Aajia dominates the cooling, is deter- 
mined by the ratio of rdiag to Aadia, which scales as / (Macc'"o) 
in self-similar disk winds for a fixed ionisation fraction. This 
scaling is smaller in our Class than in our Class 1 case by a fac- 
tor 5^ = 125, explaining the delay to heat the gas in the Class 
jet despite a 10 times lower ionisation fraction than in the Class I 
jet. 

3.1 .3. H2 survival 

The calculated H2 abundances along streamlines launched from 
1 AU and from R^ub are plotted in Fig. 6-7h, for the 3 classes 
of sources. Below, we first analyse the results in detail for the 
Class 1 case, and then discuss the Class 11 and Class cases. 

Class I case: The streamline launched just outside 
^sub - 0.2 AU suffers heavy H2 destruction (see Fig. 7h). 
Self-shielding is not operative (we assumed no molecules on 
inner, dust-free streamlines) and the FUV field is very strong. 
Hence most H2 molecules are quickly photodissociated at 



z — I AU. Balance with reformation on dust yields a final H2 
abundance ^ 10""*. 

In contrast, half of the H2 molecules are found to survive 
along the 1 AU streamline (see Fig. 6h). This may be understood 
by examining the timescales of the main formation/destruction 
processes of H2 along this streamline, plotted in Fig. 10. The 
flow timescale up to the current point, ff , is also plotted for com- 
parison and is seen to be very short (50 yrs out to 1 000 AU). 
With our assumption of a mean H2 abundance of half the lo- 
cal value, the self-shielded H2 photodissociation timescale (red 
curve in Fig. 10) is now much longer than the flow timescale. 
As a conservative check of this conclusion, we also plot in 
Fig. 10 (thick pink curve with symbols) the minimum H2 pho- 
todissociation timescale, where A^(H2) and the corresponding 
self-shielding factor are calculated assuming an H2 abundance 
on inner streamlines equal to its minimum value, as obtained on 
the streamline ro = Rsub- Even in this overly pessimistic case, 
the minimum H2 photodissociation timescale remains an order 
of magnitude longer than the flow time at all z (Fig. 10), confirm- 
ing that H2 will escape photodissociation on the 1 AU streamline 
in the Class I model. 

Only two processes are found to be fast enough to affect the 
H2 abundance on the 1 AU streamline: The first to occur is en- 
dothermic neutral-neutral reactions with O and OH to form H2O 
and H atoms (denoted collectively as O + 2H2 = H2O + 2H in 
Fig. 10): 

O + H2 ^ OH-HH-2980K (43) 
OH -H H2 ^ H2O + U- 1490 K. (44) 

The second process is collisional dissociation by H and H2. The 
timescale shortens by 6 orders of magnitude as temperature in- 
creases from 2 000 K to 4 000 K, and approaches the flow time 
around z - 2 AU, so that half of H2 is destroyed by 10 AU. 
Further out, collisional dissociation shuts off due to both the tem- 
perature decline and the density fall-off; the H2 abundance thus 
remains "frozen-in" at half its initial value (Fig. 6h). We con- 
clude that streamlines will turn from mostly atomic to mostly 
molecular around ro - 1 AU for our Class I model. 

Active Class II case: The H2 abundance shows a similar be- 
havior to the Class I case. The streamline from R^^b suffers heavy 
photodissociation, with a final abundance of ^ 2x10"^ (Fig. 7h), 
while half of the molecules survive on the 1 AU streamline 
(Fig. 6h). However, the minimum photodissociation timescale on 
the 1 AU streamline, computed with the same method as above, 
is now shorter than the flow time around z ^ 30 AU; there- 
fore, our local approximation to self-shielding might signifi- 
cantly underestimate H2 photodestruction on the Class II stream- 
line launched at 1 AU. Our H2 abundance in this particular case 
should thus be viewed conservatively as an upper limit, pending 
more detailed non-local calculations. 



Class case: Here, wind densities are so high that dust screen- 
ing against photodissociation is very effective, even on the 
streamline just outside Rsub- In addition, temperature is so low 
that oxydation reactions and collisional dissociation are slower 
than the flow crossing time, except briefly on the R^^b stream- 
line. Therefore, the Class disk wind retains 90% - 100% of 
its full initial H2 content on dusty streamlines launched from 
ro > /?sub = 0.2 AU. 
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Figure 11. Abundances relative to H nuclei of important O, C, and S-bearing species along a flow line launched at 1 AU, for our 
Class (left). Class I (middle) and Class II models (right). A color version of this figure is available in the on-line edition of this 
journal. 



3.1.4. Main O, C, and S-bearing species along tine 1 AU 
streamline 

The fractional abundances of the main O- and C-bearing neutral 
species with respect to H nuclei are plotted in Fig. 1 1 along the 
1 AU streamline for the 3 classes of sources. S and SO are also 
plotted, as they were identified recently as tracers of molecular 
jets in Class sources (Dionatos et al. 2009, 2010; Tafalla et al. 
2010). 

CO: We find that no CO destruction occurs on the Class 
streamline from 1 AU, thanks to the very efficient dust shield- 
ing provided by dense inner wind layers. In contrast, CO is 
partly destroyed by FUV photodissociation along the Class I/II 
streamlines. Once the FUV field drops down, CO reforms mainly 
through the reaction HCO^ -F e — > CO + H but the timescale is 
longer than the flow time and the final abundance per H nucleus 
is only 5 x 10"^ in the Class I jet, and 10"^ in the Class II jet. 
These latter values should be considered conservatively as upper 
limits, as CO self-shielding could be significantly overestimated 
by our approximate "local" treatment (see Section 3.2 below). 
We thus predict an abundance ratio of H2/CO > 0.5 x lO"* in all 
cases. 

H2O: When the wind temperature exceeds - 400 K, oxygen not 
locked in CO is converted to H2O through the endothermic rea- 
ctions of O and OH with H2 in eq. (43)-(44). Together with H3O+ 
recombination, this reaction is efficient enough to balance the 
water destruction processes active in Class I/II jets (FUV pho- 
todissociation, photoionisation, reaction with Hp. The asymp- 
totic abundance of water is thus similar and quite high on the 
1 AU streamline for all 3 classes\ at - 5 - 10 x 10"^. 

C, O, S atoms, OH and SO: OH and atomic C and O have low 
abundances along the Class streamline, where CO and H2O 
are well shielded. In contrast, high abundances of these species 
^ 0.5 - 1 X 10 are predicted in the outer regions of the Class I 

^ In the Class streamline, cold dust grains hold an additional H2O 
reservoir in the form of ice, of assumed abundance 10"* (I-lovvcr & 
Pineau des Forets 2003), that could be later released in the gas phase 
in shock waves. 



and Class II streamlines, following CO and H2O partial pho- 
todissociation (see Fig. 11). Concerning sulfur, it is mainly in 
the form of atomic S beyond 100 AU along the Class stream- 
line, and is mainly photoionized into S^ in the Class I/II cases, 
where the FUV field is more intense (see Fig. 8). The asymptotic 
abundance of SO is substantial in the Class only, at a level of 
3 X 10"^. These predicted characteristics for Class jets are in 
line with the relative abundance of SO to CO ^ 2 x 10"^ re- 
ported by Tafalla et al. (2010), and with the mass-flux traced by 
atomic sulfur lines being possibly as large as that inferred from 
CO (Dionatos etal. 2010). 

3.2. Effect of increasing launcli radius in the Class I model 

In Fig. 12 we illustrate in the Class I case the effect of the launch 
radius on the wind chemistry and temperature. We plot a range 
of ro of 1 - 9 AU corresponding to a range in final flow speed 
of 90 km s"' to 30 km s"'. The streamline from /?sub = 0.2 AU 
is also shown for comparison. The main effect is that the flow 
has an "onion-like" thermal-chemical structure, with stream- 
lines launched from larger radii having higher H2 abundance 
and lower temperature and ionization. The behavior with ro is 
discussed in more detail below. 

We do not present results for the interval 0.2 < ro < 1 AU in 
the Class I jet, as the minimum H2 photodissociation timescale 
there becomes similar to or less than the flow time (see discus- 
sion in Section 3.1.3). Accurate H2 abundances on these inner 
streamlines await a detailed non-local treatment of H2 shielding, 
not taken into account in the present preliminary approach. For 
the same reason, we do not present results at ro > 1 AU for the 
Class II jet, but we note that the effect of increasing ro will be 
qualitatively the same as in the Class I case (i.e. higher H2 abun- 
dance and lower temperature). We also do not present results at 
ro > 1 AU in the Class case, since most H2 and CO molecules 
survive already at a launch radius just beyond the sublimation 
radius /?sub - 0.2 AU (see Section 3.1.3). 

3.2.1 . lonisation and drift speed at various ro 

Wind streamlines are initially less ionised at larger ro 
(Fig. 12f,g), due to the larger attenuation of X-rays and FUV 
photons by intervening streamlines (Fig. 12c). At large dis- 



15 



Panoglou et al.: Molecules in protostellar MHD disk winds. I. Model 



10' 



10'' = 



10' r 



10' 




z(AU) 



lO'^ 



^ 10 



^ 10' r 



10" r 



id 

r 


1 1 i| 1 1 i| 1 11 

bulk : 


- / /7\ 
/ /•' .•' 

. ./J ,•■ ..M . . . 


drift 1 

. . .1 . . .1 . . ; 




z(AU) 




10' - 



10' 



--- 0.20 AU 

1.00 AU 

3.00 AU 

9.00 AU 

.1 . . .1 . . .1 



.1 . . .1 



z(AU) 





10 


10 




10 


12 










C\J 

X 


10 


14 










10 


16 




10 


18 



1 ' .iu-:,: 1 


"1 ' "1 ' "1 ' " 

s 












%, - 

. .1 . . .1 . . .1 . . . 




10'^ 10'^ 10° 10^ 10^ 10^ 10"* 
z(AU) 




z(AU) 



10" 
10^ 



10' 



10" - 



10"' - 



10"' - 





1 1 1 ll 1 1 ll 1 IJ 


: / / } 

— 1:1 




1 . . .1 1 . ;i . : . 


1 . . .1 . . .1 . . . 



o 
o 




z(AU) 



Figure 12. Thermo-chemical structure for a range of streamlines in our Class 1 jet; R^^b - 0.2 AU(green dashed, taken from Fig. 7), 
Class 1 jet: 0.34 AU (blue dotted), 0.58 AU (green solid), 1 AU (red solid), 3 AU (pink dotted), and 9 AU (black dashed). The arrow 
shows the effect of increasing launch radius on the plotted quantities, which are the same as in Fig. 6. CO abundances are upper 
limits, except on the ro = R^ah streamline (calculated with no self-shielding). A color version of this figure is available in the on-line 
edition of this journal. 



tances, they recombine to an asymptotic fractional ionisation of 
a few 10"^, dominated by S^. 

The drift speed does not exceed about 10 km s ' and re- 
mains less than 30% of the bulk flow speed out to ro = 9 AU 
(Fig. 12d). Therefore, magnetocentrifugal forces appear able to 
lift molecules out to disk radii of 9 AU for the typical parameters 
of Class 1 sources, in the MHD solution investigated here. 

3.2.2. Temperature and molecule survival at various ro 

The gas temperature at larger launch radii follows a similar be- 
havior to that found at ro = 1 AU (i.e. a gradual rise on a 
scale z - ro, followed by a shallow decrease), but the tem- 
perature peak is shifted to larger z (i.e. lower «h) and reaches 
a smaller peak value (see Fig. 12a). Therefore, H2 chemical 
oxydation and collisional dissociation are less efficient than for 
ro = 1 AU. As a result, the H2 abundance increases with ro and 
more than 90% of the molecules survive for launch radii > 3 AU 
(Fig. 12h). Photodissociation of H2, which was already too slow 
for ro = 1 AU, is further reduced by the additional screening 
from intervening wind streamlines and plays no role. 

In contrast, we find that photodissociation remains the ma- 
jor destruction process for CO out to ro = 9 AU. The drop in 



CO abundance towards the inner ro = R^ub streamline largely 
exceeds the factor 2 assumed in our "local" computation of the 
CO self-screening column (see Fig. 121). The plotted CO abun- 
dances for ro > Rsub should thus be considered as upper limits 
only, and a detailed non-local treatment of self-shielding will be 
needed to obtain more accurate values. 



4. Discussion 

4.1. Comparison with static disk atmospliere models 

It is instructive to compare our results with thermo-chemical cal- 
culations of irradiated static molecular disks at a radius of 1 AU. 
Glassgold et al. (2004) focussed on X-ray disk irradiation, with a 
similar X-ray flux to ours (Lx -2x 10^° erg s"', kTx - 1 keV). 
Nomura & Millar (2005) studied FUV irradiation by the star 
with standard ISM dust properties as adopted here, and a mean 
unshielded flux at a distance of 140 AU of Go - 1 000 times the 
mean interstellar flux, averaged over the interval 910-2 000 A 
(see their Fig. 4), similar to our Class 1 case (see Fig. 6f). 

We find that the thermo-chemical structure of MHD disk 
winds launched from 1 AU in Class 1 sources markedly differs 
from irradiated static disks at the same radius: H2 survives to 



16 



Panoglou et al.: Molecules in protostellar MHD disk winds. I. Model 



much greater altitudes and the gas temperature rises more grad- 
ually. We briefly discuss and explain these diff'erences below. 

4.1 .1 . H2 survival and wind sliielding 

In static disk models at 1 AU, H2 is destroyed by X-rays or FUV 
photons above z ^ 0.1 AU (Glassgold et al. 2004; Nomura & 
Millar 2005). In contrast, on the 1 AU streamline of the Class I 
disk wind, half of the gas remains molecular beyond z ^ 10 AU 
and until the end of the streamline (Fig. 6h). The survival of H2 
against photodissociation arises from several effects: 

1. The first key effect is the enhanced screening provided by 
dense inner wind streamlines compared to a hydrostatic 
flared disk geometry. As shown in Fig. 1, the self-similar 
MHD disk wind solution chosen here produces density con- 
tours that are horizontal just above the disk surface. In con- 
trast, density contours in a hydrostatic disk are strongly 
flared out, allowing a more direct illumination of the disk 
surface by the central star (see Fig. 1 of Nomma & Millar 
2005). This difference is illustrated by the large Ay values 
towards the central star in Fig. 6c; in the Class I stream- 
line at 1 AU, the "effective" dust-attenuated FUV field at 
Z - 0.3 AU is exp"^'*'' ^ 2 compared to Go - 10^ in the 
static disk of Nomura & Millar (2005) at the same (r, z). The 
peak X-ray ionisation rate is ^(Hj) - 2 x 10"" s"', instead 
of 6 X 10"^ s ' at the top of the disk atmosphere in Glassgold 
et al. (2004). 

2. The second key factor is the very short flow timescale; only 
50 years from the slow point out to 1 000 AU for the 1 AU 
streamline. Static disk models assume chemical equilibrium, 
so that H2 disappears wherever destruction reactions are 
faster than reformation. In disk winds, H2 can still survive 
this situation as long as the destruction time is longer than 
the flow time (see Fig. 10). 

4.1.2. Gas temperature profile 

In static disk models at 1 AU, FUV photoelectric heating or 
X-ray heating induce a sharp increase in gas temperature up to 
1 500 - 5 000 K within a few scale heights, at z/ro ^ 0.15 - 0.2 
(Nomura & Millar 2005; Glassgold et al. 2004). In our disk wind 
model, the rise in temperature is much more gradual, with tem- 
peratures of a few thousand K reached only around z/ro - 1. 
This difference stems both from the greater attenuation of X-rays 
and FUV photons by inner wind streamlines, and from the pow- 
erful adiabatic expansion cooling experienced by the wind after 
launching. The latter term largely exceeds FUV and X-ray heat- 
ing and is only slightly offset by drag heating (see Fig. 9), yield- 
ing a slow net heating rate. A related consequence is that tem- 
peratures above 2 000 K are reached only in wind regions of low 
density so that H2 coUisional dissociation and chemical oxyda- 
tion is limited. 

4.2. Limitations of the model and planned improvements 

The results presented here are only a preliminary attempt at ad- 
dressing the highly complex problem of the coupled ionisation, 
chemical, thermal, and dynamical structure of MHD disk winds 
from young stars. Several limitations have been introduced that 
we discuss below with our planned improvements: 

1 . We ignored feed-back of the thermal structure on the dynam- 
ics. This is in fact well justified because thermal pressure 



gradients are already negligible above the slow point where 
our integration starts, and are only dynamically important in 
the mass-loading region located deeper in the disk, around 
z - h (see Casse & Ferreira 2000). Therefore, any devia- 
tions between our computed wind temperature and the pre- 
scription adopted by Casse & Ferreira (2000) to compute the 
dynamical solution do not alter the self-consistency of the 
model. Nevertheless, we show in Appendix A that the two 
temperatures fall close to each other near the wind base. The 
existence of such "slow" warm MHD solutions might then 
be a natural outcome of the heating mechanisms present near 
the surface of an irradiated molecular disk. Thermal calcula- 
tions extending into the optically thick and resistive mass- 
loading region are under way to test this conjecture (Garcia 
et al., in preparation). 

2. We also ignored feed-back of the ion-neutral drift on the 
flow dynamics; we showed that this single fluid approxi- 
mation is well fulfilled at launch radii < 9 AU but tends 
to degrade at larger ro or small Mace, due to the scaling 
of Vin/Vp oc X(\^Y'^ M^^^. Therefore, ion-neutral drift may be 
one of the factors (together with the available magnetic flux) 
limiting the outer radial extension of molecular MHD disk 
winds in young stars. The drift will also decrease rotation 
signatures in neutral jet tracers at large ro compared to the 
single fluid prediction. Multi-fluid MHD calculations would 
be useful to further quantify this effect. 

3. We adopted a simplified treatment of shielding against pho- 
todissociation. An upcoming improvement of our model 
will be the inclusion of a non-local treatment of the self- 
screening, in order to obtain more accurate predictions of 
CO and/or H2 abundances in the Class I and Class II jets. 
Shielding against photodissociation also depends on the as- 
sumed composition of inner wind streamlines launched be- 
tween rin - 0.07 AU and /?sub- As done in Ruden et al. (1990) 
and Garcia el al. (2001a), we assumed that all wind mate- 
rial with R{0) > Rsub contributes to the Ay, i.e. that dust re- 
forms efficiently on these dense streamlines after they exit 
the sublimation surface. This may be too optimistic. On the 
other hand, when computing the minimum H2 abundance 
at ro = Rsub, we assumed that no molecules are present 
on more inner streamlines to provide self-shielding, which 
is pessimistic. Indeed, CO overtone infrared line profiles in 
Class II disks indicate the presence of rotating CO down to 
Hn - 0.04 AU (Najita et al. 2003). Another planned im- 
provement in our model will thus be to perform thermal- 
chemical modelling of dust-free streamlines launched inside 
Rsub, including H2 gas phase formation processes through 
ion chemistry and three-body reactions, and advection from 
outer radii. 

4. The dust size distribution and dust to gas ratio were assumed 
to be the same as in the ISM. This is probably adequate for 
the young Class and Class I, where the disk is young and 
still being fed by an infalling envelope, but there is evidence 
for both dust growth and dust settling starting to occur dur- 
ing the Class II phase, where the envelope has dissipated. 
This will decrease the efficiency of dust screening against 
photodissociation. On the other hand, grain growth may lead 
to a smaller dust sublimation radius, so that more molecular 
streamlines can contribute to self-shielding. Knowledge of 
the dust size stratification in turbulent magnetized disks will 
be needed to assess the net effect on H2 survival in MHD 
disk winds from evolved Class II stars. 

5. The present exploratory work considered only one particular 
steady MHD disk wind solution, with a magnetic lever arm 
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parameter A ^ 14 best reproducing current tentative rotation 
signatures in atomic jets from Class 11 (see Section 2. 1). If jet 
rotation has been overestimated, slower models with smaller 
values of A would have to be considered. These denser winds 
(higher ejection efficiency ^) would be better self-screened, 
and thus less subject to photodissociation. On the other hand, 
the effect on H2 collisional dissociation is difficult to predict 
without detailed temperature calculations: the smaller accel- 
erating Lorentz force per particle will tend to decrease ambi- 
polar heating, while the lower ionisation (due to enhanced 
screening) will tend to increase it (eq. 15). Such solutions 
are being developed and will be explored in forthcoming pa- 
pers. 

4.3. Predicted observational trends and future tests 

Our preliminary results on the thermo-chemistry of a dusty cen- 
trifugal MHD disk wind exhibit general properties that seem 
promising to explain several observed trends in the molecular 
counterparts of stellar jets. They are summarized below, and fu- 
ture observational tests are outlined. 

Concerning low-mass Class sources, the model predicts 
that dusty flow streamlines will keep most of their H2 and CO 
content at least down to the sublimation radius of ^ 0.2 AU, 
i.e. up to flow speeds of ^ 100 km s"' for M* = Q.IMq. This 
seems promising to explain the frequent detection of collimated 
fast H2 and CO in Class jets, as well as the high ratio of 
H2 to CO ^ 10"* recently estimated in one of them, HH211 
(Dionatos et al. 2010). Furthermore, tentative rotation signa- 
tures reported so far in Class jets (Lee el al. 2007, 2008) seem 
consistent with steady MHD centrifugal launching from the ex- 
pected range of disk radii ro - 0.15 - 0.6 AU, when the most 
likely poloidal speeds are adopted^. Finally, if the MHD disk 
wind operates over a radial extent from = 0.2 AU to rout 
of 1 - 9 AU, the ratio of molecular mass-flux to accretion rate 
would be ^log(rout/'"in) - 0.06 - 0.15, in the same range as ob- 
servations (Lee et al. 2007). 

As the accretion rate drops and the stellar mass increases, our 
modelling indicates that the MHD disk wind is more irradiated 
(less screening against FUV radiation) and hotter (self-similar 
scaling of Fdiag with and M^cc). Therefore, the molecular H2 
zone moves to larger launch radii, around ro > 1 AU in our 
Class 1 source, with a typical terminal speed of ^ 100 km s"' 
(cf. eq. 2). The molecular region will move even further out in 
Class II sources. This predicted trend is in line with possible 
rotation signatures reported so far in molecular Class 1/11 jets 
(HH26 and CB26), which do appear to suggest larger magneto- 
centrifugal launch radii ^3-10 AU than in Class jets, al- 
though statistics are admittedly still limited (Chrysostomou et al. 

* The launch radius of a centrifugally-driven MHD disk wind, as in- 
ferred from jet rotation signatures, scales with the poloidal speed ap- 
proximately as ro oc Vp"*'^ (Anderson et al. 2003). Lcc et al. (2007, 

2008) assumed a range Vp 100 - 200 km s"' to estimate launch radii 
in the HH211 and HH212 Class jets. The smaller Vp value seems 
more probable in HH21 1 given the very low-mass of the central source 
(0.065Mo), and yields ro = 0.15 AU. In HH212, rotation signatures 
were seen only at low radial velocities ^ 1.5-4 km s"' while SiOjet 
emission is detected up to |Ki_maxl -9-12 km s"' (Lee et al. 2008). The 
maximum radial velocities correspond very well to the proper motion 

150 km s"' of H2 knots in HH212, for the estimated jet inclination of 
^ 4° to the plane of the sky (Codella et al. 2007; Claussen et al. 1998). 
Hence, the slow SiO rotating gas probably has a lower poloidal speed 
Vp 20 - 60 km s"', yielding ro - 0.6 AU in both jet lobes (Cabrit 

2009) instead of 0.3 - 0.05 AU as estimated by Lee et al. (2008). 



2008; Launhardt et al. 2009; Cabrit 2009). This prediction also 
seems promising to explain the trend towards slow, wide mole- 
cular counterparts in older Class 11 jet sources. The striking ex- 
ample of the broad slow H2 wind encompassing the HL Tau and 
DG Tau atomic jets, and the slow CO conical flow surrounding 
the atomic jet in HH 30, both suggest a "hollow" molecular wind 
structure consistent with this idea (Takami et al. 2007; Beck et al. 
2008; Pety et al. 2006). 

A third promising trend is that gas temperature on molecular 
streamlines is predicted to increase from T ^ 700 K in Class 
jets to r ^ 2 000 - 3 000 K in Class 1/11 jets. This is in good 
agreement with temperatures ^ 300 - 1 000 K recently inferred 
in Class jets from H2 pure rotational lines studied with Spitzer 
(Dionatos et al. 2009, 2010), and ^ 2 500 K in Class l/Il jets from 
H2 ro-vibrational hnes (Takami et al. 2006; Beck et al. 2008). As 
rovibrational lines might have a greater contribution from shock- 
heated gas, H2 temperature estimates in Class 11 jets from pure 
rotational lines would be useful to confirm this result. 

While these trends are encouraging, more specific tests need 
to be carried out to validate the applicability of centrifugal 
MHD disk winds to molecular protostellar jets. Closer compar- 
ison to observations await more detailed calculations of self- 
shielding for the Class 1/11, and computations of synthetic maps 
and spectra in various tracers, taking into account beam smear- 
ing and NLTE effects (Yvart et al., in preparation). Several as- 
pects promise to be particularly discriminant: 

A first important test would be a confirmation of rotation 
in molecular jets, with a resolved spatial pattern consistent with 
centrifugal acceleration from the disk. Current measurements are 
affected by beam smearing, which can significantly lower the 
observed rotation signature compared to the actual underlying 
rotation curve (Pesenti et al. 2004). This effect could be impor- 
tant in the HH212 and HH211 Class jets in Orion-Perseus, 
whose widths < 0.2 " are smaller than the best resolution of cur- 
rent submm interferometers (Cabrit et al. 2007; Lee et al. 2007). 
The lack of resolution also makes it more difficult to recognize 
possible contamination of rotation signatures by jet precession, 
orbital motions, or shock asymmetries (see e.g. Cerqueira et al. 
2006; Correia et al. 2009 for discussions of possible artefacts). 
The ALMA interferometer will be essential to progress on this 
issue. 

A second prediction of the MHD disk wind model con- 
sidered here is that, beyond the dust sublimation radius, the 
wind base is substantially shielded from stellar photons and the 
gas heats up only gradually through ambipolar diffusion (see 
Section 3.1.2). CO rovibrational lines, excited around 1 000 K, 
would thus be formed higher up in the atmosphere (at z ^ 1 AU) 
than in static disks heated mostly by irradiation. Comparison of 
predicted CO rovibrational line profiles with observed ones in 
Class I/II sources in the near-infrared (Najita et al. 2003; Bast 
et al. 201 1) may thus provide an interesting test. The high pre- 
dicted abundance of H2O > 10"^ in warm MHD disk winds 
ejected from ^ 1 AU, at all evolutionary stages (see Fig. 11) 
is also an important characteristic that may be testable with 
Herschel or infrared observations (e.g., Kristensen et al. 2010; 
Pontoppidan et al. 2010). 

Finally, we note that the wind thermo-chemical properties 
may be locally modified by internal shock waves forming in the 
jet due to time-variability or instabilities; which molecules are 
destroyed or reformed will depend on whether the shock is of C 
("continuous") or J ("jump") type, which in turn depends cru- 
cially on the ionization fraction, magnetic field intensity, and H2 
fraction in the preshock gas (Le Bourlot et al. 2002; Flower & 
Pineau des Forets 2003). The present calculations provide for the 



18 



Panoglou et al.: Molecules in protostellar MHD disk winds. I. Model 



first time the appropriate shock "initial conditions" in an MHD 
disk wind, setting the stage to calculate the predicted thermo- 
chemistry in internal shocks, and to compare it with shock ob- 
servations. Although such developments are beyond the scope 
of the present paper, we note that SiO and CH3OH, observed to 
be enhanced by 2 to 4 orders of magnitude in Class jet knots 
(Bachiller el al. 1991; Cabrit et al. 2007; Tafalla el al. 2010), 
could in principle be released in a dusty disk wind following 
shock erosion or vaporization of grains (Gusdorf et al. 2008; 
Guilletetal. 2009, 2011). 

5. Conclusions 

We have investigated the non-equilibrium thermal-chemical 
structure of dusty streamlines in a "slow" self-similar MHD 
disk wind compatible with current observational constraints in 
atomic T Tauri microjets. We considered a range of accretion 
rates and stellar masses representative of the Class 0, Class I, and 
early Class II stages of low-mass star formation, and probing a 
decrease of 2 orders of magnitude in wind density. A detailed 
chemical network, a complete set of heating/cooling terms, and 
the effect of irradiation by stellar X-rays and FUV photons, were 
considered. Our main conclusions are the following: 

- The MHD disk wind has an "onion-like" thermal-chemical 
structure, with temperature, ionization, and radiation field all 
decreasing as the launch radius of the streamline increases. 
The wind is sufficiently ionized to accelerate neutrals out to 
disk radii of ^ 9 AU, but is sufficiently self-screened and 
cool to remain molecular on streamlines launched beyond 
some minimum radius r^. For the MHD solution explored 
here, this radius is ro - 0.2 AU (sublimation radius) for 
Class parameters, ro - 1 AU for the Class I, and ro > 1 AU 
for the Class II. 

- Key elements for the survival of H2 in MHD disk winds, as 
opposed to static irradiated disk atmospheres, are: the short 
flow crossing timescales (50 - 100 yrs on the 1 AU stream- 
line), the efficient shielding provided by inner wind stream- 
lines against stellar FUV and X-ray photons, and the strong 
adiabatic cooling which delays gas heating and limits colli- 
sional and chemical destruction. 

- Balance between ion-neutral drag and molecular line cooling 
establishes asymptotic temperatures ^ 700/2 000/3 000 K 
on the 1 AU streamline, increasing with protostellar stage. 
Temperatures at the wind slow magnetosonic point are 
within a factor 2 of the surface heating prescribed in the 
MHD solution, suggesting that the assumed high wind mass- 
loading might be naturally achievable (although thermal cal- 
culations deeper into the disk are needed to test this conjec- 
ture). 

- The chemistry of dusty irradiated MHD disk winds is 
of a complex and hybrid nature, combining the out-of- 
equilibrium "warm" chemistry of C-shocks heated by ambi- 
polar diffusion, and the photoprocesses present in PDRs and 
XDRs. Illustrations are the efficient gas-phase synthesis of 
H2O at all stages, and the substantial abundances reached by 
SH^, CH+, and H+ in the hotter and more ionized Class I/II. 
In these latter cases, OH, O, and C also reach levels similar 
to H2O while CO is strongly photodissociated. 

- The model predictions appear promising to explain sev- 
eral observational trends in molecular wind counterparts, 
namely: the presence in low-mass Class sources of 
highly collimated molecular jets with typical speed up to 
100 kms"', low temperature ^ 700 K, and estimated 



magneto-centrifugal launch radii ^ 0.2 - 0.6 AU; as well 
as the trends for higher H2 temperatures ^ 2 500 K, larger 
magneto-centrifugal launch radii, and lower collimation and 
speed as the source evolves from Class to Class I to 
Class II. 

- More discriminant tests of the applicability of MHD disk 
winds to molecular jets/winds in young stars require syn- 
thetic maps and line profiles in H2, CO, and H2O, including 
a non-local treatment of self-shielding and proper account 
for beam dilution and NLTE excitation. Such calculations 
will be presented in following papers in this series, and com- 
pared with infrared and submm observations. The effect of 
internal shocks, in particular the enhancement of species be- 
lieved to be released from grains, such as SiO and CH3OH, 
will also be explored. 

Acknowledgements. We ai'e grateful to G. Herczeg and P. Lesaft're for useful 
suggestions and discussions, and to the anonymous referee for many insight- 
ful comments that helped to improve the paper D. Panoglou, S. Cabrit, and 
P. Garcia wish to acknowledge financial and travel support through the Marie 
Curie Research Training Network JETSET (Jet Simulations, Experiments and 
Theory) under contract MRTN-CT-2004-005592. This research has made use of 
NASA's Astrophysics Data System. 



References 

Aggarwal, K. M., Benington, K. A., Burke, P. G., Kingston, A. E., & Pathak, A. 

1991, J. Ph. B, 24, 1385 
Agra-Amboage, V., Dougados, C, Cabrit, S., & Reunanen, J. 2011, A&A, 532, 

A59 

Anderson, J. M., Li, Z.-Y., Krasnopolsky, R., & Blandford, R. D. 2003, ApJ, 
590, LI 07 

Antoniucci, S., Nisini, B., Giannini, T., & Lorenzetti, D. 2008, A&A, 479, 503 

Bacciotti, R, Mundt, R., Ray, T. P, et al. 2000, ApJ, 537, L49 

Bacciotti, R, Ray, T. P, Mundt, R., Eislofi"el, J., & Solf, J. 2002, ApJ, 576, 222 

Bachiller, R. & Cernicharo, J. 1990, A&A, 239, 276 

Bachiller, R., Martin-Pintado, J., & Fuente, A. 1991, A&A, 243, L21 

Bakes, E. L. O. & Tielens, A. G. G. M. 1994, ApJ, 427, 822 

Bally, J., Reipurth, B., & Davis, C. J. 2007, Protostai-s and Planets V, 215 

Bast, J. E., Brown, J. M., Herczeg, G. J., van Dishoeck, E. F., & Pontoppidan, 

K. M. 2011, A&A, 527, Al 19+ 
Beck, T. L., McGregor, P J., Takami, M., & Pyo, T.-S. 2008, ApJ, 676, 472 
Bergin, E., Calvet, N., D'Alessio, R, & Herczeg, G. J. 2003, ApJ, 591, L159 
Blandford, R. D. & Payne, D. G. 1982, MNRAS, 199, 883 
Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., & Romanova, 

M. M. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 

479^94 

Brown, P N., Byrne, G. D., & Hindmarsh, A. C. 1989, SIAM J. Sci. Stat. 
Comput., 10, 1038 

Cabrit, S. 2007a, in Lecture Notes in Physics, Berlin Springer Verlag, ed. 

J. Ferreira, C. Dougados, & E. Whelan, Vol. 723, 21-+ 
Cabrit, S. 2007b, in lAU Symposium, Vol. 243, lAU Symposium, ed. J. Bouvier 

& I. Appenzeller, 203-214 
Cabrit, S. 2009, in Protostellar Jets in Context, ed. K. Tsinganos, T. Ray, & 

M. Stute, 247-257 
Cabrit, S., Codella, C, Gueth, F, et al. 2007, A&A, 468, L29 
Cabrit, S., Ferreira, J., & Raga, A. 1999, A&A, 343, 687 
Casse, F & Ferreira, J. 2000, A&A, 361, 1178 
Cernicharo, J. & Reipurth, B. 1996, ApJ, 460, L57+ 

Cerqueira, A. H., Velazquez, P. F., Raga, A. C, Vasconcelos, M. J., & de Colle, 

F2006, A&A, 448, 231 
Chrysostomou, A., Bacciotti, F, Nisini, B., et al. 2008, A&A, 482, 575 
Claussen, M. J., Marvel, K. B., Wootten, A., & Wilking, B. A. 1998, ApJ, 507, 

L79 

Codella, C, Cabrit, S., Gueth, F, et al. 2007, A&A, 462, L53 
Coff"ey, D., Bacciotti, F, Ray, T. P, Eisloffel, J., & Woitas, J. 2007, ApJ, 663, 
350 

Coff"ey, D., Bacciotti, F, Woitas, J., Ray, T. P, & Eislolfel, J. 2004, ApJ, 604, 
758 

Combet, C. & Ferreira, J. 2008, A&A, 479, 481 

Correia, S., Zinnecker, H., lUdgway, S. T., & McCaughrean, M. J. 2009, A&A, 
505, 673 

D'Alessio, P, Calvet, N., Hartmann, L., Lizano, S., & Canto, J. 1999, ApJ, 527, 
893 



19 



Panoglou et al.: Molecules in protostellar MHD disk winds. I. Model 



Dalgamo, A., Yan, M.. & Liu, W. 1999, ApJS, 125, 237 
Davis, C. J., Hodapp, K. W., & Desroches, L. 2001a, A&A, 377, 285 
Davis, C. J., Nisini, B., Takami, M., et al. 2006, ApJ, 639, 969 
Davis, C. J., Ray, T. P., Desroches. L., & Aspin, C. 2001b. MNRAS, 326, 524 
Davis, C. J., Stem, L., Ray, T. P., & Chrysostomou, A. 2002, A&A, 382, 1021 
Davis, C. J., Whelan, E., Ray, T. P, & Chrysostomou, A. 2003, A&A, 397, 693 
Dionatos, O.. Nisini, B., Cabrit, S., Kristensen, L., & Pineau Des Forets, G. 2010, 
A&A, 521, A7+ 

Dionatos, O., Nisini, B., Garcia Lopez, R., et al. 2009, ApJ, 692, 1 

Draine, B. T. 1978, ApJS, 36, 595 

Draine, B. T. & Bertoldi, F. 1996, ApJ, 468, 269 

Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89 

Draine, B. T. & Malhotra, S. 1993, ApJ, 414, 632 

Ferreira, J. 1997, A&A, 319, 340 

Ferreira, J., Dougados, C, & Cabrit, S. 2006. A&A, 453, 785 
Flower, D. R., Le Bourlot. J., Pineau des Forets, G.. & Cabrit, S. 2003, MNRAS. 
341,70 

Flower, D. R. & Pineau des Forets, G. 2003, MNRAS, 343, 390 
Flower, D. R., Pineau des Forets, G., & Hartquist, T. W. 1985, MNRAS, 216, 
775 

Flower, D. R., Pineau des Forets, G., & Walmsley, C. M. 2007, A&A, 474, 923 
Garcia, P J. V., Cabrit, S., Ferreira, J., & Binette. L. 2001b. A&A, 377, 609 
(Paper II) 

Garcia, P J. V., Ferreira, J., Cabrit, S., & Binette, L. 2001a, A&A, 377, 589 
(Paper I) 

Giannini, T.. McCoey, C, Caratti o Garatti, A., et al. 2004, A&A, 419, 999 
Glassgold, A. E., Mamon, G. A., & Muggins, R J. 1989, ApJ, 336, L29 
Glassgold, A. E., Mamon, G. A., & Muggins, R J. 1991, ApJ, 373, 254 
Glassgold, A. E., Najita, J., & Igea. J. 1997, ApJ, 480, 344 
Glassgold, A. E., Najita, J.. & Igea, J. 2004, ApJ, 615, 972 
Gredel, R., Lepp, S., & Dalgarno, A. 1987, ApJ, 323, L137 
Gredel. R.. Lepp, S., Dalgarno, A., & Herbst, E. 1989, ApJ, 347, 289 
Gudel, M., Telleschi, A.. Audard, M., et al. 2007, A&A, 468, 515 
Gueth, F & Guilloteau, S. 1999, A&A, 343, 571 

Guillet, v., Jones, A. P, & Pineau Des Forets, G. 2009, A&A, 497, 145 
Guillet, v., Pineau Des Forets, G., & Jones, A. P 201 1, A&A, 527, A123+ 
Guilloteau, S., Bachiller, R., Fuente, A., & Lucas, R. 1992, A&A, 265, L49 
Gullbring, E., Calvet, N., Muzerolle, J., & Martmann, L. 2000, ApJ, 544, 927 
Gusdorf, A., Pineau Des Forets, G., Cabrit, S., & Flower, D. R. 2008, A&A, 490, 
695 

Hartigan, P, Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736 

Herczeg, G. J., Linsky, J. L., Valenti, J. A., Johns-Krull, C. M., & Wood, B. E. 

2002, ApJ, 572,310 
Herczeg, G. J.. Linsky, J. L., Walter, F M., Gahm, G. F, & Johns-Krull, C. M. 

2006. ApJS, 165, 256 
Hirano, N., Liu, S.-Y. Shang, M., et al. 2006, ApJ, 636, L141 
Hollenbach. D. & McKee, C. F 1979, ApJS, 41, 555 
Hummer, D. G. 1963, MNRAS, 125, 461 

Imanishi, K., Nakajima, M., Tsujimoto, M., Koyama, K., & Tsuboi, Y. 2003, 
PASJ, 55, 653 

Johns-Krull, C. M., Valenti, J. A., & Linsky, J. L. 2000, ApJ, 539, 815 
Konigl, A. & Pudritz, R. E. 2000, Protostars and Planets IV, 759 
Kristensen, L. E.. Visser, R.. van Dishoeck, E. F, et al. 2010, A&A, 521, L30+ 
Laor, A. & Draine, B. T. 1993, ApJ. 402, 441 

Launhardt, R., Pavlyuchenkov, Y, Gueth, F, et al. 2009, A&A, 494, 147 
Lavalley-Fouquet, C. Cabrit, S., & Dougados, C. 2000, A&A, 356, L41 
Le Bourlot, J., Pineau des Forets, G., & Flower, D. R. 1999, MNRAS, 305, 802 
Le Bourlot, J., Pineau des Forets, G.. Flower. D. R., & Cabrit, S. 2002, MNRAS, 
332, 985 

Lee, C.-R, Mo, R, Beuther, H., et al. 2006, ApJ, 639, 292 

Lee, C.-R, Mo, R, Hirano, N., et al. 2007, ApJ, 659, 499 

Lee, C.-R, Ho, R T. P, Bourke, T. L., et al. 2008, ApJ, 685, 1026 

Lee, C.-R, Ho, R T. P, Palau, A., et al. 2007, ApJ, 670, 1 188 

Lee, C.-R, Stone, J. M., Ostriker, E. C, & Mundy, L. G. 2001, ApJ, 557, 429 

Lee, H.-H., Herbst, E., Pineau des Rorets, G., Roueff, E., & Le Bourlot, J. 1996, 

A&A, 311,690 
Lizano, S., Heiles, C, Rodriguez, L. F, et al. 1988, ApJ, 763 
Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977. ApJ, 217, 425 
McCaughrean, M. J., Rayner, J. T., & Zinnecker, H. 1994, ApJ, 436, L189 
Najita, J., Carr, J. S.. & Mathieu. R. D. 2003, ApJ, 589, 931 
Neufeld, D. A. & Kaufman, M. J. 1993. ApJ, 418, 263 
Nisini, B., Bacciotti, F, Giannini, T.. et al. 2005, A&A, 441, 159 
Nisini, B., Codella, C, Giannini, T., et al. 2007, A&A, 462, 163 
Nomura, H. & Millar. T. 2005, A&A, 438, 923 
Pesenti, N., Dougados, C, Cabrit, S., et al. 2004, A&A, 416, L9 
Pety, J., Gueth, F. Guilloteau, S.. & Dutrey, A. 2006, A&A, 458, 841 
Pinto, C. & Gain, D. 2008a, A&A, 484 
Pinto, C. & Galh, D. 2008b, A&A, 492 



Podio, L., Bacciotti, F, Nisini, B., et al. 2006, A&A, 456, 189 
Podio, L., Medves, S., Bacciotti, R, Eisloffel, J., & Ray, T. 2009, A&A, 506, 779 
Pontoppidan, K. M., Salyk, C, Blake, G. A., & Kaufl, H. U. 2010, ApJ, 722, 
L173 

Pyo, T.-S., Hayashi, M., Kobayashi, N., et al. 2006, ApJ, 649, 836 
Pyo, T.-S., Kobayashi, N., Hayashi, M., et al. 2003, ApJ, 590, 340 
Raga, A. & Cabrit, S. 1993, A&A, 278, 267 

Raga, A. E., Williams, D. A., & Lim, A. J. 2005, RevMexAA, 41, 137 
Rapp, D. & Englander-Golden, R 1965, J. Chem. Ph., 43, 1464 
Rawlings, J. M. C, Drew, J. E., & Barlow, M. J. 1993, MNRAS, 265, 968 
Rawlings, J. M. C, Williams, D. A., & Canto, J. 1988, MNRAS, 230, 695 
Ray, T., Dougados, C, Bacciotti, R, Eisloflfel, J., & Chrysostomou, A. 2007, in 

Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 231-244 
Roberge, W. G., Jones, D., Lepp, S., & Dalgamo, A. 1991, ApJS, 77, 287 
Rodriguez, L. F, Escalante, V, Lizano, S., Canto, J., & Mirabel, I. R. 1990, ApJ, 

365, 261 

Ruden, S. R, Glassgold, A. E., & Shu, R H. 1990, ApJ, 361, 546 
Safier, R N. 1993, ApJ, 408, 1 15 

Santiago-Garcfa, J., Tafalla, M., Johnstone, D., & Bachiller, R. 2009, A&A, 495, 
169 

Schunk, R. W. 1975, Planet. Space Sci., 23, 437 
Shang, H., Allen, A., Li, Z.-Y, et al. 2006, ApJ, 649, 845 
Shang, H., Glassgold, A. E., Shu, F H., & Lizano, S. 2002, ApJ, 564, 853 
Shang, H., Shu, F H., & Glassgold, A. E. 1998, ApJl, 493, L91 
Shu, R H., Lizano, S., Galh, D., Cai, M. J., & Mohanty, S. 2008, ApJ, 682, L121 
Shu, F H., Najita, J. R., Shang, H., & Li, Z.-Y. 2000, Protostars and Planets IV, 
789 

Shu, S., Gain, D., Lizano, S., & Cai, M. 2007, in lAU Symposium, Vol. 243, 

lAU Symposium, ed. J. Bouvier & I. Appenzeller, 249 
Stabler, S. W. 1988, ApJ, 332, 804 
Suzuki, T. K. & Inutsuka, S.-i. 2009, ApJ, 691, L49 

Tafalla, M., Santiago-Garci'a, J., Hacar, A., & Bachiller, R. 2010, A&A, 522, 
A91-I- 

Takami, M., Beck, T. L., Pyo, T.-S., McGregor, R, & Davis, C. 2007, ApJ, 670, 
L33 

Takami, M., Chiysostomou, A., Ray, T. R, et al. 2004, A&A, 416, 213 

Takami, M., Chiysostomou, A., Ray, T. R, et al. 2006, ApJ, 641, 357 

Taylor, S. D. & Raga, A. C. 1995, A&A, 296, 823 

Tielens, A. G. G. M. & Hollenbach, D. 1985, ApJ, 291, 722 

van Dishoeck, E. F. 1988, in Rate Coefficients in Astrochemistry, ed. T. Millar 

& D. Williams (Dordrecht, Boston: Kluwer Academic Publishers), 49 
Wolk, S. J., Hamden. Jr. R R., Placcomio, E.. et al. 2005, ApJS, 160, 423 



Appendix A: Comparison with the assumed heating 
in the dynamical solution 

The MHD solution adopted in the present work has been ob- 
tained with a prescribed self-similar enthalpy injection func- 
tion peaking near the disk surface, representing a tiny fraction 
/ = 8 X 10""* of the accretion power (Casse & Ferreira 2000). 
This small surface heating is sufficient to obtain denser, slower 
disk winds with mass fluxes and poloidal and rotation speeds 
consistent with observations (see Section 2.1), but the nature of 
the heating process is not yet identified. 

In this context, it is interesting to compare the thermal pres- 
sure assumed in the underlying MHD accretion-ejection solution 
with that obtained from our thermal-chemical calculations where 
the various known heating, cooling and chemical processes are 
taken into account. Since the gas density p is prescribed by the 
solution, the relevant quantity to compare thermal pressures is 
just the ratio T Ifi (with ^ being the mean molecular weight). 
This comparison is shown in Figure A. 1 along a 1 AU stream- 
line for our Class 0, Class I, and Class II models. 

The value of Tdyn/yUdyn in the dynamical solution first un- 
dergoes an increase from the disc surface at z = /i = 0.03ro 
to about Zt - lO/i and then an adiabatic decrease. Although the 
initial rise is not exactly followed by our calculations, it is inter- 
esting that the difference between the two curves is less than a 
factor 2 near the slow magnetosonic point (our starting integra- 
tion point for the thermo-chemistry). This suggests that the warm 



20 



Panoglou et al.: Molecules in protostellar MHD disk winds. I. Model 




10' 



-I I— J I 1— Ll I 1— L 



1^7— rrr I 



10"^ 10"'' 10° 10^ 10^ 10^ 
z/rn 




10' 



_l l—U I L-Ll I LJJ I I— U I L_L 



10"^ 10"'' 10° 10^ 10^ 10^ 
z/fn 



pressure will not affect the wind dynamics. Although the thermal 
pressure gradient did play an important role in the mass loading 
(done within the resistive MHD disk layers), it plays no role in 
the jet acceleration and collimation itself which is mainly due to 
magnetic forces. Thermal pressure gradients will eventually start 
affecting the jet dynamics after the jet refocusses towards the 
axis. For this reason we have not extended our thermal-chemical 
calculations beyond the recolhmation point of the MHD solu- 
tion. 



Figure A.l. Comparison of T/jU (in K) as calculated from the 
heating/cooling terms in our thermal-chemical models (red solid 
curves) with the same ratio as assumed in the dynamical MHD 
solutions of Casse & Ferreira (2000), in dashed green, for the 
1 AU streamUne and 3 stars of differing accretion properties. 
Note that in the dynamical solution, Tdyn is calculated assuming 
a fully ionised gas, with jUdyn = 0.5. A color version of this figure 
is available in the on-line edition of this journal. 



gas surface temperatures required to load dense, slow MHD disk 
winds compatible with observational constraints might be natu- 
rally achievable under a wide range of accretion rates, especially 
as MHD "turbulent/wave" heating will take place in the deeper 
resistive disk layers, and wind mass-loading could be further en- 
hanced via upward motions induced by MRI turbulence (vSuzuki 
& Inutsuka 2009). 

The discrepancy between computed and assumed T can 
exceed a factor 10 at z, = lO/i, and a factor 100 beyond. This 
discrepancy has however no relevant impact on the dynamical 
behavior of the solution. Indeed, at z,, the ratio of plasma pres- 
sure to magnetic pressure in the dynamical solution is already 
down to only a few 10"^, so even a large discrepancy in thermal 
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